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ABSTRACT 

Increasingly large samples of galaxies are now being discovered at redshifts z ~ 5-6 and higher. 
Many of these objects are inferred to be young, low in mass, and relatively unreddened, but detailed 
analysis of their high quality spectra will not be possible until the advent of future facilities. In this 
paper we shed light on the physical conditions in a plausibly similar low mass galaxy by presenting 
the analysis of the rest-frame optical and UV spectra of Q2343-BX418, an L* galaxy at z = 2.3 
with a very low mass-to-light ratio and unusual properties: BX418 is young (< 100 Myr), low mass 
(M* ~ 10 9 Mq), low in metallicity (Z ~ 1/6 Z Q ), and unreddened (E(B -V)~ 0.02, UV continuum 
slope (3 — —2.1). We infer a metallicity 12 + log(0/H) = 7.9±0.2 from the rest-frame optical emission 
lines. We also determine the metallicity via the direct, electron temperature method, using the ratio 
O III] AA1661, 1666/[0 III] A5007 to determine the electron temperature and finding 12 + log(0/H) = 
7.8 ±0.1. These measurements place BX418 among the most metal-poor galaxies observed in emission 
at high redshift. The rest- frame UV spectrum, which represents ~ 12 hours of integration with the 
Keck telescope, contains strong emission from Lya (with rest-frame equivalent width 54 A), He II 
A1640 (both stellar and nebular), C III] AA1907, 1909 and O III] AA1661, 1666. The C IV/C III] ratio 
indicates that the source of ionization is unlikely to be an AGN. Analysis of the He II, O III] and C III] 
line strengths indicates a very high ionization parameter logt/ ~ — 1, while Lya and the interstellar 
absorption lines indicate that outflowing gas is highly ionized over a wide range of velocities. It remains 
to be determined how many of BX418's unique spectral features are due to its global properties, such 
as low metallicity and dust extinction, and how many are indicative of a short-lived phase in the early 
evolution of an otherwise normal star-forming galaxy. 

Subject headings: galaxies: abundances — galaxies: evolution — galaxies: high-redshift 



1. INTRODUCTION 

Many of the key questions in the study of galaxy for- 
mation and evolution are also the most difficult to ad- 
dress, because they center on the properties of the very 
faint galaxies in the early universe. Significant numbers 
of such ga l axies are now being confirmed at z ~ 5-6 (e.g . 
Yan et all [20051 IMcLure et al.l [200l lAndo et all |2007; 
Vanzella et al.ll2009HStark et al.ll2009ll2010[ ). and the re- 
cent installation of the Wide Field Camera 3 (WFC3) 
on the Hubble Space Telescope has resulted in a flurry of 
candidate ga laxies at 7 < z < 10 dBouwens et al.ll2010a : 
lOesch et all [2010l : lYan et al.l [20091 : iBunker et al.l 12009 : 
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IMcLure et al.l 120101 ) . The general properties of these 
galaxies can be inferred from their broadband spectral 
energy distributions, but as a result of their faint contin- 
uum magnitudes, spectral information — when available 
at all — is limited to strong Lya emission or low signal- 
to-noise composite spectra. Because of this flux limit, 
significant progress in the understanding of the detailed 
properties of these galaxies is unlikely to occur before 
the era of JWST and the extremely large ground-based 
telescopes. 

An alternative and immediately accessible approach is 
to identify galaxies at somewhat lower redshifts that are 
likely to be similar to these high redshift objects. Galax- 
ies at z ~ 2 are still bright enough for detailed study with 
current technology, and have the further advantage that 
the full set of rest-frame optical emission lines are observ- 
able in the near-IR from the ground. If we wish, then, to 
identify z ~ 2 galaxies that are most likely to illuminate 
the properties of galaxies in the very early universe, what 
features should we look for? Intuitively, we expect that 
such galaxies will be young, low in mass and metallic- 
ity, and relatively unreddened. Recent results bear this 
out. At increasing redshifts, the UV-continuum slope 
(3 (fx tx A' 3 ) of galaxies becomes bluer; the most lumi- 
nous g alaxies have ~ —1 -5 at z ~ 2, and ~ —2.4 at 
z ~ 6 (jBouwens et alj |2009). At a given redshift, fainter 
galaxies also show bluer slopes. At z ~ 7-8, galaxies 
are found to have ~ —2 to —3 (jBouwens et alj[2010bt 
iFinkelstein et a l. 2009). Given the limited data available 
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on these galaxies, the origin of such blue slopes is diffi- 
cult to determine. iBouwens et al.l (|2010bl ) suggest that 
the bluest slopes of the faintest galaxies cannot be pro- 
duced by normal stellar population s at nonzero metallici- 
ties, while IFinkelstein et al.l (|2009( ) argue that recourse to 
Population III stars, a top-heavy initial mass function, or 
other unusual mechanisms is not necessary to account for 
the blue color. Instead, they find that the galaxies' ob- 
served SEDs are best explained by a young (< 100 Myr), 
low metallicity (Z < 0.1 Z at 1 a confidence) stellar 
population with little to no dust extinction and stellar 
masses of 10 8 -10 9 M©. We are unlikely to find galaxies 
dominated by Populat ion III stars at z ~ 2, b ut the set 
of properties found by IFinkelstein et al.l (2009) may well 
be approximated by some unusual lower redshift objects. 

In this paper we present the analysis of the rest-frame 
ultraviolet and optical spectra of one such highly un- 
usual galaxy. Q2343-BX418 is an L* galaxy at z = 
2.3 se lected by the stan dard rest-frame UV color cri- 
teria ()Steidel et al.l 120041) . but it is nevertheless among 
the youngest and lowest stellar mass continuum-selected 
z ~ 2 galaxies identified to date, with rest-frame B-band 
mass-to-light ratio M/L = 0.03. More uniquely, BX418 
also shows extremely low dust obscuration and low 
metallicity. We use the UV and optical line ratios in com- 
bination with photoionization modeling to show that the 
properties of this galaxy are significantly different from 
those of more typical z ~ 2 galaxies, and from local star- 
burst galaxies as well. Our study relies on an extremely 
deep UV spectrum: given spectra of sufficient quality, ob- 
servations in the rest-frame UV are sensitive to dust con- 
tent, age, galaxy-scale outflows, metallicity, ionization 
parameter, and the initial mass function. In order to ob- 
tain the relatively high signal-to-noise necessary for this 
analysis, detailed studies of the UV spectra of high red- 
shift galaxies have so far been restri cted primarily to the 
composite spectra of many galaxies (|Shaplev et al.ll2003l) 
or to galaxies whose luminosities ha ve fortuitously been 
boost e d by gravitational l e nsing dPettini et all 
2002t iCabanac et all 120081: iQuider et all 120091 
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Dessauges-Zavadskv et al.ll2010n . The result is that these 



studies have primarily served to illuminate the properties 
of typical galaxies at z ~ 2-3. Here we take advantage 
of ~ 15 hours of integration with the Keck telescope to 
study an object that is considerably younger, bluer, and 
lower in mass than these more typical galaxies. 

The paper is organized as follows. We describe the ob- 
servations and data reduction in Section [21 the construc- 
tion of a composite spectrum for comparison in Section 
O and the global properties of BX418 in Section 0j_in- 
cluding the stellar population properties (Section 14. ip . 
the star formation rate (Section I4.2j) . and the kinemat- 
ics and gas mass (Section I4.3j) . In Section [5] we discuss 
the suite of photoionization models used to match the 
observed line ratios. We focus on the rest-frame optical 
spectrum in Section [5J including constraints on metallic- 
ity and ionization parameter. The rest-frame UV spec- 
trum is the subject of Section [7j Section 17.11 deals with 
the interstellar absorption and lyct emission, and Sec- 
tion 17.21 describes the He II emission. In Section 17.31 we 
discuss nebular emission from carbon and oxygen and 
resulting constraints on the C/O abundance ratio. The 
stellar wind lines are described in Section 17.41 and the 
weak fine structure emission lines in Section 17.51 We 



discuss all of these results in Section [51 and summarize 
our conclusions in Section [9l We use a cosmology with 
H = 70 km s" 1 Mpc" 1 , fl m = 0.3, and n A = 0.7, 
and assume a IChabrierl ([20031 ) initial mass function and 
a solar oxygen abun dance of 12 + log(0/H) = 8.69 
([Asplund et al.H2009D . 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Optical Spectroscopy 

The z = 2.3048 galaxy Q2343-BX418 (a = 
23:46:18.571, 6 = +12:47:47.379; J2000) was observed 
in the course of th e z ~ 2 galaxy survey described by 
iSteidel et all ([20041 ). In September and October 2003, 
we created four slitmasks in the Q2343 field for deep ob- 
servation wi th the low Res olution Imaging Spectrom- 
eter (IRIS; lOke et all 119951 ) on the Keck I telescope. 
BX418 appeared on two of these masks, and was ob- 
served for 5.5 hrs on 24-26 September 2003 and for 6.1 
hrs on 26-28 October 2003, for a total of 11.6 hrs of in- 
te gration. The observa tions were conducted as described 
bv lSteidel et all (|2004l ). We used the 400 groove mm" 1 
grism blazed at 3400 A for the IRIS blue channel, and 
divided the beam with a dichroic at 6800 A. In the red 
channel we used the 400 groove mm^ 1 grating blazed at 
8500 A, though in the case of BX418 all of the scientif- 
ically useful information falls in the blue channel. Each 
mask was observed at the parallactic angle, for approxi- 
mately two hours each night. The deep masks contained 
a total of 77 objects, receiving integrations ranging from 
5.5 to 22.8 hrs. Inspection of the spectra showed many to 
be of high quality but simi lar to the average UV spectra 
of galaxies at z ~ 2-3 (e.g. iShaplev et al.ll2003l) ; BX418, 
however, contains several distinctive spectral features. 

Data reduction wa s done accord i ng to standard pro- 
cedures described bv ISteidel et all ((2003). The spectra 
from each of the two masks containing BX418 were re- 
duced and extracted independently, and then these ID 
spectra were combined to produce the final spectrum. 
Flux calibration was checked by integrating the spec- 
trum over the G bandpass using the IRAF task "sbands" 
and comparing with the broadband magnitude. The two 
measurements differed by only 0.03 mag, so no correc- 
tion to the flux calibration was applied. The data re- 
duction procedures also included the production of a 2D 
variance spectrum for each slit on each mask observed. 
These were scaled by the number of exposures of each 
mask, extracted to ID using the same apertures used for 
the science spectra, and combined using standard error 
propagation to construct a la error spectrum which was 
used to determine uncertainties. The spectrum of BX418 
covers wavelengths ~ 3400-6400 A, or ~ 1000-1930 A in 
the rest frame. The slit width was 1.2", but because the 
galaxy is compact and the seeing was good, the spectral 
resolution is primarily determined by the angular size of 
the galaxy, which was ~ 0.85". The resulting resolution 
is ~ 400 km s _1 , and ranges from ~ 550 km s _1 at the 
blue end of the spectrum to ~ 300 km s _1 at the red 
end. 

2.2. Near-IR Spectroscopy 

BX418 wa s also observed with NIRSPEC 

(|McLean et all I1998D on the Keck II telescope, in 
order to obtain infrared spectra of the rest-frame optical 
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emission lines. The if- band obse r vation of Ha was 
included in the survey of lErb et al.l (|2006d ); these data 
were obtained on 13 June 2004, and consisted of five 15 
minute integrations using the NIRSPEC-6 filter. We 
also observed BX418 in the H (for [O III] AA5007, 4959 
and H/3) and J (for [O II] A3727) bands in September 
2004. The iJ-band observations were conducted on 
6 September 2004 using the NIRSPEC-5 filter, and 
consisted of four 15 minute integrations. Similarly, we 
obtained four 15 minute J-band integrations using the 
NIRSPEC-3 filter on 7 September 2004. As measured 
from the sky lines, the spectral resolution is 15 A (~ 200 
km s -1 ) in the K band, 10 A (~ 180 km s -1 ) in H, 
and 10 A (~ 240 km s" 1 ) in J. The data were reduced 
using the standard procedures described by lErb et all 

poap . 

When comparing observations taken at different times 
and with different spectral configurations, relative (and 
absolute) flux calibration is always a concern. As dis- 
cussed bv lErb et al] (|2006d ). NIRSPEC flux calibration 
is subject to significant uncertainties, primarily because 
of slit losses. In the case of BX418, however, we have rea- 
son to believe that these uncertainties are minimal. As 
described below, Ha emission from BX418 has also been 
observed with OSIRIS, the near-IR integral field spectro- 
graph with laser guide star adaptive optics on the Keck II 
telescope. The two flux measurements agree within 3%, 
well below the uncertainties on both observations. The 
OSIRIS observations also show that the nebular emission 
arises from a very small region with radius 0.8 kpc, or 
0.1" at z — 2.3. Thus the spatial extent of the galaxy in 
the NIRSPEC observations is determined almost entirely 
by the seeing, which was ~ 0.5" for all observations. The 
slit width was 0.76", and we therefore expect slit losses 
to be minimal. Conditions were photometric for all ob- 
servations. As noted above, the flux calibration of the 
rest-frame UV spectrum was checked by comparison with 
the broadband G magnitude and found to agree; we are 
therefore confident in comparisons between the LRIS and 
NIRSPEC spectra as well. 

2.3. Optical and IR Imaging 

We also make use of deep optical, near-IR and mid-IR 
images. The Q2343 field was imaged in the U n , G and 
R s bands with the Large Format Camera (LFC) on the 
Pal omar 200-inch telesc ope in August 2001, as described 
by ISteidel et all (pOOl . Near-IR J and K band ob- 
servations were conducted with the Wide-Field Infrared 
Camera (WIRC. IWilson et al.ll2003ft . also on th e Palomar 
telesco pe. These observations are described bv lErb et al.l 
(|2006d ). Finally, the Q2 343 field was obse rved at 4.5 and 
8.0 /xm with the IRAC (jFazio et al.ll2004h instrument on 
the Spitzer Space Telescope, and at 24 /urn with MIPS. 
We use only the 4.5 /im IRAC data, as BX418 was not 
detected at 8 /im. The 4.5 /im image had an integration 
time of 18000 sec and a 3a depth of ~0 16 nJy. The 
IRAC image was reduced as described bv lShaplev et al.l 
(2005). BX418 was also undetected in the MIPS 24 /tin 
image, which had an effective integration time of 9.6 hrs 
and a 3<r depth of ~15 /tJy. The MIP S data were reduced 
as described bv lReddv etail (j2006b[ ). 

A portion of the Q2343 field including BX418 was ob- 
served using the F160W filter (A ~ 1.54/tm, tracing rest- 



frame 4700 A at the redshift of BX418) of the WFC3/IR 
camera on board the Hubble Space Telescope (HST) on 
UT 15 June 2010, as part of a survey of star-forming 
galaxies in deep extragalactic fields (GO-11694). The in- 
tegration time was 9x900 sec, for a total of 8100 sec. 
The details of the HST observations and data reduc- 
tion are described in full by Law et al. (in preparation); 
in brief, the HST images were drizzled to a pixel scale 
of 0.08 arcsec pixel -1 with a point spread function of 
FWHM = 0.19" (corresponding to 1.6 kpc at the red- 
shift of BX418), and reach a 5a limiting magnitude of 
27.9 AB in a 0.4" diameter aperture. 

3. COMPARISON SAMPLE 

In order to compare the rest-frame UV spectrum of 
BX418 with that of more typical galaxies at z ~ 2, 
we have constructed a composite spectrum of all of the 
galaxies in the z ~ 2 sample for which we have deter- 
mined stellar population parameters, and which have 
spectra of high enough signal-to-noise for reliable red- 
shift determination. The composite is the average of 
1186 spectra of 966 galaxies (because some galaxies were 
observed more than once). It was constructed in the 
same fashion as the z ~ 3 Lyman break galaxy com- 
posite spectrum described bv lShaplev et alj (12003). and 
is very similar to that spectrum in appearance. The 
galaxies included in the z ~ 2 composite have the fol- 
lowing medians and semi-interquartile ranges: redshift 
z mcd = 2.2±0.3, stellar mass M*, med = 9.7±8.9xl0 9 M , 
and age 571±538 Myr, assuming constant star formation. 
The median reddening is E(B-V) me d = 0.17±0.06. The 
properties of BX418 are given in Section [4~T1 for compar- 
ison. 
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Figure 1. The G — TZ vs. 1Z — Kab color-color diagram, in 
which BX418 is compared with 1116 spectroscopically confirmed 
z ~ 2 galaxies selected by the BX, BM and MD color criteria 
ISteidel et a l. 2003, 20Ql). 



4. GLOBAL PROPERTIES 

Broadband magnitudes for BX418 are given in Ta- 
ble [1] from U n to 4.5 /im. Note that the spectral en- 
ergy distribution is nearly flat, with AB magnitudes 
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Table 1 

Photometry 





G 


n 




^AB a 


^AB b 


4.5 fim 


24.3±0.1 


23.9±0.1 


24.0±0.1 


24.5±0.4 


23.04±0.04 


23.7±0.4 


23.8±0.5 



Note. — All magnitudes are AB. 
a HST F160W magnitude, with no correction for [O III] AA4959, 5007 and H/3 
emission. The corrected magnitude is 23.49±0.08. 
b Uncorrected for Ha emission, which accounts for 0.4 ±0.1 mag. 



mAB — 24 ±0.5 in all observed bands after correction for 
line emission. These colors are highly unusual, as shown 
in Figure [TJ Here we present the G — TZ vs. TZ — K color- 
color diagram for 1116 spectroscopically confirmed z ~ 2 
alaxies selected by the BX, BM and MD color criteria 
Steidel et all [2001 [200l . BX418 is shown by the large 
circle in the lower left corner of the diagram. There are 
few objects in this corner, and BX418 lies at the extreme 
such that galaxies bluer in G — TZ are redder in 7Z — K, 
while galaxies bluer in 1Z — K are redder in G — 1Z. In 
physical terms, the blue G — TZ color is indicative of low 
extinction, but not necessarily a young stellar popula- 
tion, while the blue 1Z—K color implies a low stellar mass 
and a young age as well, given the high star formation 
rates of these galaxies. Young, low mass, and unreddened 
galaxies can therefore be expected to fall in the lower left 
corner of the color-color diagram; they are few in number 
because young galaxies usually have redder UV slopes 
UShaplev et al.ll200ll [2005 IReddv et al.ll2010h . From the 
G — TZ color, we find a UV continuum slope = —2.1, 
placing BX418 at the extreme blue end of the distri- 
butio n of luminous z ~ 2 galaxies (jReddv et al.1 l2006bL 
2010). Note that at z = 2.3 Lya falls between the U n 
and G bands, so the observed colors are not affected by 
the strong Lya emission. However, the H and K magni- 
tudes are affected by the strong rest-frame optical emis- 
sion lines; combined [O III] AA4959, 5007 and H/3 emission 
contributes 0.45 mag to the observed .ff-band, while Ha 
emission contributes 0.4 mag to the if-band magnitude. 

BX418 is undetected in deep 8 and 24 /an imaging, 
with 3<7 upper limits of tuab > 22.5 at 8 /im and 
/24 [im < 34 /ijy. Most AGN in the UV-selected z ~ 2 
galaxy sample have a rising, pow er law SEP at 8 an d 
24 /im, with / 24 ^m > 100 /xJy (|Reddv et all l2006aj ): 
BX418's nondetection in these bands therefore argues 
against the presence of an obscured AGN. 

In Figure [5] we present a high resolution rest frame op- 
tical image of BX418, taken with the Wide Field Camera 
3 on the Hubble Space Telescope, using the F160W filter. 
This demonstrates that the galaxy is very compact, with 
a half light radius of 1.5 kpc; for comparison, the average 
half light radius of galaxies in the UV-selected sample is 
~ 3 kpc. BX418 shows a slight elongation alo ng the N-S 
axis (also seen in the OSIRIS observations of lLaw et al.l 
2009! described in Section [43] below), but otherwise has 
no distinctive morphological features. 

4.1. Stellar Population Modeling 

Stellar population properties are determined by mod- 
eling of th e broadband SEP in the usual fashion, em- 
ploying the lBruzual fc Charlotl ()2003l ) population synthe- 
sis models and fo ll owing met hods previously de scribed 
by IShaplev et all (pOOl and lErb et al.l (|2006d ). We 
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Figure 2. Rest-frame optical image of BX418, taken with the 
Wide Field Camera 3 on the Hubble Space Telescope. The area 
shown is 3" X 3" , and the grey circle at lower right shows the 
FWHM of the 0.19" PSF. BX418 is compact, with a half light 
radius of 1.5 kpc. 



quantify the uncertainties in the modeling using a se- 
ries of 10,000 Monte Carlo simulations, as described by 
IShaplev et al.l (|2005l ). These simulations perturb the 
galaxy's photometry in accordance with the photomet- 
ric errors, and compute the best fitting model for the 
perturbed photometry, including a variety of star for- 
mation histories. We find that 91% of the 10,000 trials 
favor a young, low mass galaxy with age < 150 Myr; 
the remaining 9% allow a more massive and significantly 
older stellar population. Based on the strong spectral 
evidence that the galaxy is young and metal poor, we 
adopt the younger and less massive solution, taking as 
the stellar mass, age, and E(B — V) the mean values 
of these quantities in the 91% of trials which give this 
general solution. Uncertainties are then given by the 
standard deviation. With the exception of a somewhat 
younger age (38 Myr vs. 100 Myr), these values are not 
significantly different from the best fit constant star for- 
mation model given by the observed photometry. We 
determine the stellar population parameters in this way 
in order to obtain uncertainties which are not skewed by 
the tail of older solutions disallowed by the spectrum. 
We have also checked these results using a newer set of 
population synthesis models which offer an improved de- 
scription of thermally pulsating AGB stars (Chariot & 
Bruzual 2007; S. Chari ot, private communic ation), with 
the SMC rather than iCalzetti et ail ((2000) extinction 
law, and with lower metallicity models of 0.4 Zq rather 
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than solar. None of these variations significantly change 
the results. We find M * = 9 x 10 8 M , an age of 38 Myr, 
and E(B — V) = 0.02. These values and their associated 
uncertainties are also given in Table [2] 

In addition we have assessed the probability of a hidden 
older stellar population by fitting t wo component ste llar 
population models, as described bv lErb et al.l (|2006cD . A 
general two component model, in which the age and mass 
of both components are allowed to vary arbitrarily, pro- 
vides a good fit to the data but increases the total mass 
by at most a factor of two. A maximally massive model, 
in which the optical and IR light is fit by a maximally 
old stellar population and the UV residuals are then fit 
with an episode of current star formation, increases the 
total stellar mass by a factor of ten. However, this model 
is a significantly poorer fit to the observed photometry, 
and overpredicts the star formation rate by a factor of 
> 2 relative to that predicted by independent observa- 
tions of Ha with both NIRSPEC (long slit spectroscopy) 
and OSIRIS (integral field spectroscopy). We therefore 
consider a significant underlying older stellar population 
to be unlikely. 

4.2. Star Formation Rates and Timescales 

We use the NIRSPEC Ha luminosity to determine a 
star formation rate of 15 ± 2 M Q yr _1 , where no ex- 
tinction correction has been applied but the uncertain- 
ties include both statistical flux errors and a range of 
extinction < E(B — V) < 0.05. This is a factor of 
~ 2 lower than the average Ha-determined s tar forma- 
tion rate of z ~ 2 galaxies (Er b et al.l 1200 6b) . We also 
calculate the SFR from the 1500 A luminosity, finding 
SFRuv = 9 Mg yr _1 , again with no extinction correc- 
tion app lied. For b o th calc ulations we use the prescrip- 
tions of iKennicuttl (|1998af ). It is well known that Ha 
emission and the UV continuum luminosity probe differ- 
ent timescales; Ha emission arises in H II regions around 
the most massive stars, and therefore traces the nearly in- 
stantaneous star formation rate, while the UV luminosity 
continues to rise for the first ~ 100 Myr. In a very young 
galaxy, we would therefore expect Ha to give a higher 
star formation rate than the UV continuum. This is in- 
deed the case in BX 418, where SFRuv/S FR Hq = 0.65. 
Using a Starburst99 (jLeitherer et all 1999) constant star 
formation model matched as closely as possible to the 
properties of BX418, we predict the UV continuum lu- 
minosity as a function of time, and find that it reaches 
~ 65% of its final value at an age of ~ 15 Myr. This is 
lower than our estimated age, but within the uncertain- 
ties. It is also comparable to the estimated dynamical 
time of ~ 25 Myr (see below). Applying our best fit 
E(B - V) = 0.02 to both the Ha and UV fluxes, we find 
SFRuv/SFRh q — 0.8, which occurs at an age of ~ 40 
Myr; this is well-matched to our best fit age. However, 
we caution that while all of these numbers are consistent 
with the young age we derive for BX418, uncertainties 
in the extinction mean that the ratio of SFRs alone is 
insufficient to reliably constrain the age to < 100 Myr. 
Applying a correction of E(B — V) = 0.05, still within 
the uncertainties, to both the Ha and UV fluxes results 
in SFRs that are equal within the photometric and spec- 
troscopic uncertainties. 

We can also place a lower limit on the SFR from the 



Lya emission, from which we find SFRL yQ = 6 M Q yr" 1 . 
A comparison of the SFRs derived from Ha and Lya sug- 
gests a Lya escape fraction of ~ 40%, as also discussed 
in Section [7X21 

4.3. Kinematics, Gas Mass and Gas Fraction 

The Ha emission from BX418 has also been observed 
with the integral field spectrograph OSIRIS on the Keck 
II telescope using laser guid e star adap t ive op tics, as part 
of the survey described by lLaw et aLl (|2009f ) . The mea- 
sured Ha flux was in excellent agreement with that de- 
rived from our NIRSPEC observations. The OSIRIS ob- 
servations indicate that the Ha emission is very compact, 
covering 2.2 ±0.6 kpc 2 with effective radius r = 0.8 ±0.1 
kpc. BX418 also shows a low mean velocity dispersion 
Cmcan = 61±17 km s _1 (where cr moa n is the flux- weighted 
mean of the velocity dispersions of individual spatial pix- 
els) and a slight velocity gradient given by « s hear = 23±12 
km s _1 over a distance of ~ 1 kpc; u s hear is defined as half 
of the maximum difference in velocity between any two 
postions along a pseudo-slit oriented to maximize the ve- 
locity gradient. Thus like most low to average mass z ~ 2 
galaxies studied to date, BX418 appears to be kinemati- 
cally dominated by random motions rather than rotation 
(jLaw et all [2009t iForster Schreiber et all 120091) . From 
the Ha kinematics and size, lLaw et aLl (|2009l ) derive a 
dynamical mass Md yn = 3 x 10 9 M Q . We also use the 
size and velocity dispersion to estimate the dynamical 
time, tdyn ~ 2r/cr ~ 25 Myr. This is again compara- 
ble to the age derived from SED fitting, and provides an 
approximate lower limit to the age of the galaxy. 

The flux and spatial extent of the Ha emission ob- 
served by OSIRIS can be used to estimate the mass 
of cold gas associated with star formation. We use 
the inversion of the K ennicutt- Schmidt star form ation 
law (|Kennicuttlll998bD given by lErb et all (|2006cD . and 
find M gas = 2.5 x 10 9 M Q . We combine this with the 
stellar mass M± — 9 x 10 s M Q to find baryonic mass 
Mbar = 3.4 x 10 9 Mq, in good agreement with the dynam- 
ical mass, and gas fraction fi = M gas / (M+ + M gas ) = 0.7. 
This high gas fraction is consistent with the galaxy's 
young age and low metallicity, and is comparable to the 
high g as fractions infer red for other young galaxies at 
z ~ 2 (|Erb et al J 120063) . 

5. PHOTOIONIZATION MODELING 

In the following sections we compare the rest-frame 
optical and UV emission line ratios with a suite of pho- 
toionization models in order to obtain constraints on 
metallicity, the ionization parameter and ionization cor- 
rection factors, and the relative abundances of C and 
O. The models were constr ucted with ve r sion 8.00 of 
Cloudy, last described by iFerland et ail (jl998f ). We 
employ an input ion izing spectrum from Starburst99 
(jLeitherer et al.l H999). assuming continuous star forma- 
tion at 15 M Q yr" 1 , an age of 50 Myr, and a metallicity 
of Z = 0.004, or ~ 0.3 Z Q , the closest available metal- 
licity to that measured for BX418. We use the Paul- 
drach/Hillier stell ar atmospheres, w hich have ionizing 
fluxes discussed bv lSmith et al.l (|2002l ). Uncertainties as- 
sociated with these ionizing fluxes are discussed in more 
detail in Section [7721 We consider metallicities of 0.05, 
0.1, 0.2, 0.3 and 0.5 Z Q . We use Cloudy's default H II 
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Table 2 

Stellar Population, Gas and Kinematic Properties 





Age a 


E(B - V) a 


SFR Ha b 


M gas c 


M iyn d 


(1O 9 M ) 


(Myr) 




(M Q yr" 1 ) 


(10 9 M Q ) 


(10 9 M Q ) 


0.9 ±0.4 


38 ±23 


0.02 ± 0.03 


15 ±2 


2.5 


3 



a Stellar mass, age and E(B — V) from SED modeling as described in the 
text. 

b Star formation rate from He?. No extinction correction has been ap- 
plied, but uncertainties include both random flux error and the effects of 
extinction up to E(B - V) < 0.05. 
c Gas mass inferre d from the Kennicutt-Schmidt law. 
d ILaw et all I pOOSh 



region abundance set, and in addition reduce the abun- 
dances of C and N by factors of 0.85, 0.65, 0.45, 0.25, 
0.15, and 0.05 for each of the metallicities adopted. The 
default C/O and N/O ratios are 0.75 and 0.17 respec- 
tively; the resulting reduced C/O ratios are 0.65, 0.49, 
0.34, 0.19, 0.11 and 0.04, and the reduced N/O ratios are 
0.15, 0.11, 0.079, 0.044, 0.026 and 0.0088. 

The ionization parameter represents the ratio of ioniz- 
ing photon to particle density, and is defined in dimcn- 
sionless form as 



U = 



Q(H) 



(1) 



where Q(H) is the number of ionizing photons emitted 
per unit time, r Q is the assumed distance from the ion- 
ization source to the cloud, and n(H) is the total hydro- 
gen density. Because neither the particle density nor the 
ionization parameter can be measured independently for 
BX418, we instead construct a set of models with rea- 
sonable values of n(H) and r Q that result in ionization 
parameters ranging from \ogU = —3.3 to \ogU = —0.2. 
We examine the variation of the line ratios of interest 
as a function of ionization parameter, metallicity, and 
C/O and N/O abundance, and choose the models best 
matching all of the observations. As described in detail 
in the following sections, the properties of BX418 can 
only be reproduced by models with low metallicity, a 
significantly reduced C/O ratio, and an extremely high 
ionization parameter. 

6. THE REST-FRAME OPTICAL SPECTRUM 

The J, H and K band observations described in Sec- 
tion [2] cover the approximate rest frame wavelength 
ranges 3450-4200, 4550-5400, and 5900-7000 A, respec- 
tively, and therefore include the lines [O II] A3727, H/3, 
[O III] AA4959, 5007, Ha, and [N II] A6583. Upper 
limits on [O II] and [N II] and fluxes of the remaining 
lines are reported in Table [31 and the spectra are shown 
in Figure El We use the strongest of these lines, Ha 
and [O III] A5007, to determine the systemic redshift of 
BX418, z sys = 2.3048. In this section we use the ratios 
of all of these lines to place constraints on metallicity, 
ionization parameter, and extinction. 

6.1. Metallicity and Ionization Parameter 

The gas phase metallicity of BX418 can be deter- 
mined from the ratios of the rest-frame optical emission 
lines, using the various familiar strong-line metallicity 



diagnostics^ We begin by obtaining an upper limit 
on the metallicity using the ratios of [N II] /Ha (N2) 
and ([N II] /Ha) /(\Q III] /H/3) (03N2) as calibrated by 
IPettini fc Pageil (l200l . We obtain 12 + log(0/H) < 8.1 
from N2, and 12 + log(0/H) < 8.0 from 03N2. These 
upper limits place BX418 on the lower branch of the 
double- valued R23 me tallicity calibration, as does the 
[N II] /Ha ratio itself ([Kewlev fc Ellison! l200l . Given 
this constraint, we can then use the R23 method to de- 
termine metallicity. This is the ratio of the oxygen lines 
[O III] AA4959, 5007 and [O II] A3727 to H/3, and is 
perhaps the most commonly used stro ng line me tallicity 
diagnostic. We use the calibr ation o f McGaugh (1991), 
in the analytical form given bv lKobulnickv et al.l JT999). 
This form incorporates the dependence of R23 on the 
ionization parameter via the inclusion of the ionization 
parameter-sensitive ratio [O III]/[0 II]. Our lower limit 
log([0 III]/[0 II] )> 1.1 indicates a high ionization pa- 
rameter of at least log U ~ — 2. In determining the metal- 
licity, we consider the range 1.1 < log([0 III]/[0 II]) < 2, 
where the upper limit comes from the maximum value 
of this ratio in our Cloudy photoionization models, cor- 
responding to a low metallicity Z = 0.05 Zq and a very 
high ionization parameter log U ~ — 1. We do not ap- 
ply an extinction correction, which would be negligible 
given the very low reddening implied by the SED mod- 
eling. Incorporating both the uncertainty in the [O II] 
flux and random error in the fluxes of the other lines, we 
find 12 + log(0/H) = 7.9 ± 0.2, or Z = 1/6 Z Q . This 
places BX418 among the most metal-poor star-forming 
galaxies known at high redshift (much lower abundances 
are found among damped Lya absorp tion systems, which 
show a wider range i n met allicity; e.g. IPettini et alJll997l : 
iProchaska & W olfe 199§). This is shown in Figure IH 
which places BX418 on the mass-metallicity relation at 
z ~ 2. 



In Section [7731 we also derive the metallicity via the di- 
rect, electron temperature method, using the ratio O III] 
AA1661, 1666/[0 III] A5007 to determine the electron 
temperature. We find 12 + log(0/H) = 7.8 ± 0.1, in ex- 
cellent agreement with the strong line metallicity derived 
above. From later work on the UV spectrum (Section 
7.2p . we also find that reproducing all of the line ratios 
requires a very high ionization parameter of log U ~ — 1. 
We discuss this extreme ionization parameter in Section 



10 Unfortunately the T e -sensitive auroral line [O III]A4363 falls 
outside of our observational windows, though given the faintness of 
this line we would be unlikely to detect it even in the low metallicity 
case. 
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Figure 3. Rest-frame optical spectra in the K, H and J bands, taken as three separate observations with NIRSPEC. The thin blue lines 
show the la error spectra, and the vertical dashed lines mark the locations of the strong emission lines. From left to right in the K band: 
[N II] A6548 (undetected), Ha, [N II] A6583 (undetected), and [S II] AA6716, 6731 (undetected); in the H band, H,3, [O III] A4959, and 
[O III] A5007; and in the J band, [O II] A3727 (undetected). 



6.2. Extinction and the Ha/H/3 Ratio 

We have a further estimate of dust extinction in BX418 
in the Ha/H/3 ratio. In principle this provides a more ro- 
bust measurement than the SED fitting, but in practice 
we are limited by the low S/N of the H/3 line. For case 
B recombination, Ha/H/3 = 2.8; from our measured Ha 
and H/3 fluxes, we find Ha/H/3 = 3.1 ± 0.4. This is con- 
sistent with zero extinction and with the low value of 
E(B — V) estimated from the SED fitting, but allows for 
greater reddening as well. We conclude that the S/N of 
the Ha/H/3 ratio is too low to provide additional con- 
straints on the amount of dust extinction. 

6.3. The [N II] /Ra vs. [O III]/RP Diagnostic Diagram 
and the [N II]/Ha Ratio 



The [N II] /Ha vs. [O III] /H/3 line ratio diagram is one 
of several diagnostic tools commonly used to gain insight 
into the physical conditions in galaxies and distinguish 
betw e en star-forming galaxies and AGN (Baldwi n et al.l 
[19811: IVeilleux k Osterbrocl [198711 . In this section we 
use the upper limit on [N II] and our measurements of 
the other three lines to place BX418 on this diagram, and 
our suite of photoionization models to interpret the re- 
sults. The diagram is shown in Figure [SJ where the large 
black circle shows the limit on the position of BX418. 
For context, the small grey points show ~ 97, 000 objects 
from the Sloan Digital Sky Survey; the left branch shows 
the location of star- forming galaxies, while AGN extend 
upward and to the right. The dotted and dashed lines 
show empirical and theoretical divisions be tween galaxies 
and AGN ffrom lKauffmann et al.ll2003l and lKewlev et all 
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Table 3 

Emission Lines 



Ion 


\ L a 

^rcst 


W b 


Flux 


FWHM C 


Av d 


Instrument 




(A) 


(A) 


(10 — 17 erg s _1 cm -2 ) 


(km s _1 ) 


(km s" 1 ) 




Lya 
Si II* 


1215.6701 


54. 0±1. 2 


29.3±0.4 


840±17 


307±3 


LRIS-B 


1264.738 


0.4±0.1 


0.2±0.06 




220±72 


LRIS-B 


Si II* 


1309.276 


0.5±0.2 


0.2±0.07 




206±71 


LRIS-B 


C TT* 
O 11 


looO. / Uo 


n 74-n 1 
U. / ±v. 1 


U.oztU.UD 




iy / ztoi 


f DTC D 


Si II* 


1533.431 


< 0.4 


< 0.2 






LRIS-B 


He II 


1640.42 


2.7±0.2 


0.8±0.1 


612±64 


-17±55 


LRIS-B 


O m] 


1660.809 


1.0±0.2 


0.3±0.07 


< 276 


193±57 


LRIS-B 


O in] 


1666.150 


1.3±0.2 


0.4±0.08 


235±93 


-2±41 


LRIS-B 


C m] 


1906.683, 1908.734 d 


7.1±0.4 


1.4±0.1 


225±36 


81±21 


LRIS-B 


[O ii] 


3727.09, 3728.79 




< 1.9 






NIRSPEC 


H/3 


4862.721 


44±9 


2.6±0.3 


< 102 


-46±19 


NIRSPEC 


[O m] 


4960.295 


93±13 


5.4±0.3 


120±23 


-47±10 


NIRSPEC 


[O m] 


5008.239 


285±26 


16.7±0.2 


161±7 


-9±4 


NIRSPEC 


Ha 


6564.614 


450±247 


8.0±0.2 


155±11 




NIRSPEC 


[N ii] 


6585.27 




< 0.3 






NIRSPEC 



a Vacuum rest wavelength 

b Rest frame equivalent width, computed from the spectrum for the LRIS-B observations and from the 
comparison of the line flux and broadband magnitude for the rest-frame optical emission lines. No broadband 
magnitude is available for the other lines observed with NIRSPEC. Where the lines are undetected, 3cr limits 
are given. 

c Full width at half maximum, corrected for instrumental resolution. If the uncorrected line width minus its 
error is less than the instrumental resolution, we report a lcr upper limit of the corrected FWHM plus its 1<t 
error. We do not report widths for the weak UV fine structure lines, which are unresolved. 
d Velocity offset relative to Ha emission. Velocity differences of H/3 and [O III] A4959 with respect to [O III] 
A5007 are due to sky line residuals. 

c We assume a wavelength of 1907.709 A for the blended C III] doublet. 




10 11 
Log (M„„/M s ) 

Figure 4. BX418 and the mass-metallicity relation for z ~ 2 
galaxies. All measuremen t s are either from the [N II] /Ha cali- 
bration of Pcttini & Pagcl (2004 ) or have been converte d to this 
method using the prescriptions of[Kcwlcy & Ellison (2008). BX418 
is indicated by the large magenta circle at lower left. Other 
data points are from Law et al. 2009 (L09), Forstcr Schrcibcr et al. 
1006] (FS06) , IGenzel et al.H200a . 120081 (G06, 08). I Yuan fc Kewlevl 
2009 (Y09), IStark et al.1120081 (S08). and the composite spectra of 
Erb ct al. 2006a (E06). The horizontal dotted line shows the solar 
oxygen abundance. 



120011 respectively) . The ope n red stars show ot her galax- 
ies at z ~ 2 as discussed bv lErb et al.l (|2006al ). 

Next we plot the line ratios predicted by the Cloudy 
models, using filled circles color-coded by ionization pa- 
rameter. We restrict the models shown to the set with 
Z = 0.2 Z Q for clarity; for higher metallicities, the curves 
move up and to the right. We see that the [O III] /H/3 



ratio increases with increasing ionization parameter (be- 
cause the fraction of oxygen in the +2 ionization state 
increases with ionization parameter, while the strength 
of H/3 depends on the star formation rate). [N II] /Ha, on 
the other hand, decreases with increasing ionization pa- 
rameter (as the N + /N +2 ratio decreases with increasing 
ionization parameter). Unsurprisingly, the [N II] /Ha ra- 
tio also depends on the nitrogen to oxygen abundance 
ratio N/O. While we have no constraints on N/O in 
BX418, it is expected that the nitrogen abundance will 
be depressed relative to oxygen in young galaxies, be- 
cause N is produced by intermediate mass stars as well 
as by the massive stars that produce O and is there- 
fore released into the ISM more slowly. This effect has 
been seen in the lensed z = 2.1 galaxy MS1512-cB58, in 
which N is underabundant with respect to O by a factor 
of ~ 3 (|Pettini et al.l [20021 ) . a s well as in local ga laxies 
and in damped Lya systems (|Pettini et al.| [2008). We 
therefore plot curves for several values of N/O in Figure 
El The rightmost curve uses the default N/O ratio of 
Cloudy's H II region abundance set; proceeding to the 
left, e ach curve further re duces the abundance of N (for 
cB58. iPettini et al.l pOOl find log (N/O) = -1.89). 

We use this diagram to make several points. First, the 
[N II] flux is likely to be well below our lower limit; we are 
unlikely to detect it even with significantly deeper data, 
especially if the nitrogen abundance is low with respect 
to oxygen, as is likely. We also note that under these 
extreme conditions of very high ionization parameter and 
(possibly) low N/O ratio, the N2 metallicity calibration 
would provide a significant underestimate of the oxygen 
abundance. This is not unexpected, as the local sample 
from which the relation is calibrated contains no galaxies 
with such an extreme ionization parameter P^l 

11 The R23 diagnostic we use to determine the metallicity in Sec- 
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Figure 5. The [N Il]/Ha vs. [O IIl]/H/3 diagnostic diagram. BX418 is indicated by the large black circle, and the small grey points 
show ~96,000 objects from the Sloan Digital Sky Survey. Colored points are predictions from photoionization models, coded by ionization 
parameter as labeled at lower left. The ionization parameter increases upward and to the left along the curves. All models shown have 
20% of the solar oxygen abundance, and curves are shown for six different values of N/O, illustrating the effect of decreasing nitrogen 
abundance on the [N Il]/Hct ratio. For higher overall oxygen abundances, the curves move up and to the right. The open red stars show 
other galaxies at z ~ 2 from Erb ct al. (2006a]); only objects with de tections of all four lines are s hown. The dotted an d dashed lines show 
empirical and theoretical divisions between galaxies and AGN, from Kauffmann ct al. (2003) and Kcwlcy ct al. (2001) respectively. 



This very low [N II] /Ha ratio places BX418 far from 
the region of the diagram occupied by AGN, suggesting 
that significant AGN contribution to the ionizing flux of 
BX418 is unlikely. However, most of the AGN shown 
in Figur e [5] are considerably higher in metallicity than 
BX418. I Groves et all (|2006l ) have studied emission line 
diagnostics in lower metallicity AGN, finding that the 
[N II] /Ha ratio is much more metallicity-dependent than 
[O III]/H/3, and that therefore low metallicity AGN move 
to the left on the diagram, with relatively little vertical 
shift. Given BX418's low metallicity and probable low 
nitrogen abundance, its position on the diagnostic dia- 
gram alone is not sufficient to rule out the presence of 
an AGN. However, as we have argued elsewhere, the lack 

tion l6.ll does not suffer from the same issue, as it incorporates the 
ionization parameter; we have used the best-fitting Cloudy models 
to verify that the R23 index reproduces the input metallicity of the 
models within our uncertainties. 



of strong C IV emission indicates a radiation field softer 
than that seen in AGN (Section[7]), and the nondetection 
of BX418 in deep 8 and 24 /im images argues against the 
presence of an obscured AGN. 

Next we observe that the few z ~ 2 galaxies for which 
all four emission lines have been detected lie significantly 
to the right of the probable location of BX418 (and from 
the s eque nce of local star-forming g alaxies; see lLiu et al.1 
120081 and lBrinchmann et al.ll2008bl ). This is true by defi- 
nition, given our limited sensitivity to [N II] , but stacking 
of spectra shows that most sta r- forming z ~ 2 galax- 
ies have log ([N II] /Ha) > -1 (|Erb et all l2006al ) . This 
much higher [N II] /Ha ratio confirms that average z ~ 2 
galaxies (while still having ionization parameters higher 
than are seen in typical star-forming galaxies in the local 
universe) have a lower ionization parameter and higher 
metallic ity than BX418, in agreement with ot her obser- 
vations (|Hainline et al.ll2009t lErb et alj|2006al) . 
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7. THE REST-FRAME UV SPECTRUM 

Rest-frame UV spectra of galaxies contain a myriad of 
features sensitive to dust content, age, galaxy-scale out- 
flows, metallicity, i onization parameter , and the initial 
mass function (e.sz. iShaplev etaTI [2001 . BX418 is not 
unusual among z ~ 2 galaxies in showing strong Lya 
emission and strong absorption from high ionization in- 
terstellar lines, but it is rare in that the low ionization 
absorption lines are weak or undetected, and in its ex- 
tremely strong emission from He II A1640, C III] AA1907, 
1909 and O III] AA1661, 1666. The spectrum is shown in 
full in Figure |6l and with an expanded vertical scale for 
detail in Figure [7J 

Before proceeding with detailed discussion of the spec- 
trum, we ask whether these strong emission lines indicate 
the presence of an AGN. We compare BX418's emission 
line ratios (see Table \S§ with those of the co mposite spec- 
trum of 16 narrow-line AGN presented by iSteidel et al.l 
(2002). These AGN have a mean redshift (z) = 2.67, 
and were selected i n the course of th e z ~ 3 Lyman 
break galaxy survey (jSteidel et al.ll2003h . The composite 
AGN spectrum shows many of the same emission lines 
as BX418, but is notably different in its much stronger 
C IV emission. The hard AGN ionizing spectrum is ev- 
ident in the ratio C IV/C III] ~ 2; in contrast, BX418 
has C IV/C III] ~ 0.3, indicative of a significantly softer 
radiation field. The composite AGN spectrum also has 
relatively strong C IV in comparison with Lya, with 
C IV/Lya ~ 0.25, while BX418 has C IV/Lya - 0.01. 
Given the weakness of the C IV emission, which can be 
fully accounted for by nebular and P Cygni stellar emis- 
sion (see Section 17. 4[) , significant AGN contribution to 
BX418's ionizing flux is unlikely. 

7.1. Outflowing Gas 

The strongest features in the UV spectrum of a typi- 
cal high redshift galaxy — a variety of interstellar absorp- 
tion lines, and Lya emission and/or absorption — reflect 
the kinematics of po werful, galaxy-scale outflows of gas 
(jShapley et al.ll2~003h . In nearly all such spectra, the in- 
terstellar absorption lines are blueshifted with respect to 
the systemic redshift of the galaxy, while Lya emission 
(if present) is redshifted. The standard model invoked to 
explain this pattern is of a nearly spherical outflow; the 
blueshifted absorption lines arise from the gas approach- 
ing the observer, on the near side of the galaxy, while 
the resonantly scattered Lya photons are most likely to 
escape the galaxy in the direction of the observer when 
they are shifted away from their resonant frequency by 
scattering off the receding gas on the far side of the 
galaxy, resulting in a redshift. The interstellar absorp- 
tion lines therefore provide a reasonably straightforward 
probe of the kinematics of the outflowing gas, subject 
to the uncertainties imposed by limited spectral resolu- 
tion and partial covering fraction. Lya emission also re- 
flects the kinematics of the outflow, but the complexities 
of Lya radiative transfer make its interpretation much 
more difficult. Several recent models explore the connec- 
tions between Lya (and int erstellar) line profiles and the 
kinematics of the gas (e.g. iVerhamme et al.l l2006. 2008; 
ISteidel et al1l20Toh . 

7.1.1. Interstellar Absorption 



A detailed comparison of the UV spectrum of BX418 
and the composite spectrum of ~ 1000 z ~ 2 galaxies is 
shown in Figure [7J The composite spectrum is charac- 
terized by strong low-ionization absorption lines from Si, 
O, C, Fe and Al, originating in outflowing neutral and 
low ionization gas as described above. In BX418 these 
lines are weak or absent. We adopt a 3cr threshold for 
detection of these lines; applying this threshold, the only 
low ionization line with a statistically significant detec- 
tion is Si II A1527. Other lines that appear to be weakly 
present in Figure[7J such as Si II A1260, O I + Si II A1303 
and C II A1334, fall just below this threshold. 

In typical galaxies at z ~ 2-3, the low ionization lines 
are strongly saturated, indicating that the equivalent 
width is determined by the velocity range and cover- 
ing fraction of the absorbing gas. Although the lines 
are barely present in BX418, we can nevertheless ob- 
tain interesting constraints by comparing the two Si II 
lines at 1260 and 1527 A. On the linear, optically thin 
part of the curve of growth, the equivalent width of the 
line W oc NX 2 f, where N is the column density, A is 
the wavelength of the transition, and / is the oscillator 
strength. For the two Si II lines, the oscillator strengths 
of the blue and red tr ansitions are 1.18 and 0.133 respec- 
tively ([Morton! 120031 ). so that in the optically thin case 
W(1260)/W(1527) = 6.0. We measure W = -0.7 A for 
Si II A1527, while the 1260 line is undetected with a 3ct 
limit Wq > —0.5 A. In the optically thin case we would 
find VT(1260) ~ —4 A, strongly ruled out by the data. 
We also use the Cloudy models to assess the strength 
of any predicted emission from the Si II lines, finding 
that such emission is at most ~ 10% of the strength of 
the observed absorption, too weak to be detected in the 
spectrum or affect the measured line ratios. We therefore 
conclude that in spite of their weakness, the low ioniza- 
tion lines are saturated, as in other galaxies at z ~ 2-3. 

The spectrum also contains high ionization absorp- 
tion lines from Si IV and C IV. In contrast to the low 
ionization lines, these are as strong or stronger than 
their counterparts in the composite spectrum. As with 
the Si II lines discussed above, the Si IV doublet at 
AA1393, 1402 A can be used to obtain constraints on 
the optical depth of the absorbing gas. The oscillator 
strengths of the members of the doublet (see TableS]) im- 
ply W(1393)/Vy(1402) = 2.0 in the optically thin case; 
the line strengths in the composite spectrum are con- 
sistent with this rati o, as they are i n the z ~ 3 LBG 
composite studied bv IShaplev et "all (|2003h . In BX418, 
however, the two Si IV lines are of approximately equal 
strength, indicating that the high ionization lines are sat- 
urated. Although the composite spectra show that this 
is not the usual case, BX418 is not the only such exam- 
ple; the high ionization lines are also saturated in sev- 
eral of the lensed galaxies as well, t hough these objects 
also have strong low ionization lines dPettini et alJl2000b 
lOuider et aHl200ftlDeIsauges-Zavadskv et al.ll2010D . We 
have also checked the predicted strength of nebular emis- 
sion from the two Si IV transitions: this emission is ex- 
pected to be stronger than the Si II emission considered 
above, but for the range of ionization parameters BX418 
requires (log U > — 2 to log U ~ —1), emission from 
the two members of the doublet is approximately equal, 
leaving the absorption line ratio unaffected. 
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Figure 6. The rest-frame UV spectrum of Q2343-BX418, shown in full to emphasize the strong Lya emission. The strongest features are 
marked with vertical dotted lines: from left, Lya emission, the interstellar high ionization absorption lines Si IV A1394, Si IV A1402, and 
C IV A1549, He II A1640 emission, O III] A1663 emission, and C III] A1909 emission. The spectrum is not normalized; note the flatness of 
the UV slope. Overplotted for comparison (thin red line) is the composite spectrum of 966 galaxies with a mean redshift (z) = 2.2. 



The equivalent width of a strongly saturated line is 
determined by the velocity range and covering fraction 
of the absorbing gas; if the line is resolved, the transi- 
tion describes the covering fraction of g function 
of velocity. Given the fairly low resolution and limited 
S/N of the BX418 spectrum, it is difficult to distinguish 
between these two effects, and detailed modeling of the 
kinematics of the outflowing gas is not feasible. How- 
ever, the available information suggests that the range 
of outflow velocities is typical of z ~ 2-3 galaxies. The 
Si IV lines have a mean blueshift of ~ 150 km s _1 , and 
show asymmetric absorption extending to approximately 
—700 km s _1 relative to the systemic velocity defined by 
the nebular emission lines (Figure [8j left two panels); 
both the centroid and range are cha racteristic of inter- 
stellar absorption in z ~ 2 galaxies ( Pcttini et alj[2000t 
IQuider et alJ2009ll20iarSteidel et al.ll2010D . The S/N of 



the low ionization lines is too low to make a reliable mea- 
surement of the velocity range, but the mean blueshift of 
the Si II A1527 line is similar to that of the Si IV lines, 
and the low and high ionization lines typically show sim- 
ilar kinematics. The Lya emission line (discussed below) 
also shows a very wide range of velocities. This would 
suggest that covering fraction is the primary factor driv- 
ing the relative strengths of the low and high ionization 
lines, and that the covering fraction of gas producing the 
high ionization lines is significantly higher than that re- 
sponsible for the low ionization lines. Given the extreme 
ionization parameter (see Section lT.2[) . this result is not 
surprising. 

7.1.2. Lya Emission 

BX418 is a strong Lya emitter, with rest-frame equiv- 
alent width Whya = 54 A. The Lya line is also unusu- 
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Figure 7. The normalized UV spectrum shown on an expanded vertical scale. We overplot the normalized z ~ 2 composite in red for 
comparison. Lines in both spectra are identified. The features marked "int" at 3859 and 4926 A are intervening Lyo and C IV absorp tion 
from the z = 2.176 foreground galaxy Q2343-BX442, with an impact parameter of 140 kpc. This system is discussed in Section 17.1.31 
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emission. The two left panels show the Si IV AA1393, 1402 lines, and the right panel shows Lya. 



A YOUNG, UNREDDENED, LOW METALLICITY GALAXY AT HIGH REDSHIFT 



Table 4 

Interstellar Absorption Lines 



Ion 


\ a 
^rcst 


f h 


W C 


Av d 




(A) 




(A) 


(km s- 1 ) 


Si II 


1260.4221 


1.18 


> -0.5 




O 1° 


1302.1685 


0.04887 


> -0.3 




Si ll e 


1304.3702 


0.094 


> -0.3 




C II 


1334.5323 


0.1278 


> -0.5 




Si iv 


1393.7602 


0.5140 


-1.1 ± 0.1 


-153 


Si iv 


1402.7729 


0.2553 


-1.3 ±0.1 


-161 


Si ii 


1526.7070 


0.133 


-0.7 ± 0.1 


-91 


C iv f 


1548.204 


0.1908 


-5.1 ± 0.2 


-583 


C iv* 


1550.781 


0.09522 


-5.1 ±0.2 


-583 


Fe ii 


1608.4509 


0.058 


> -0.6 




Al H« 


1670.7886 


1.833 


> -0.7 





a Vacuum rest wavelength. 

b Oscillator strength, from Morton (2003). 

c Rest frame equivalent width. Where the lines are undetected, 
3<r limits are given. 

d Velocity offset relative to Ha emission. 
c Blended. 

' Blended; Wo is given for the blend, and Av assumes that the 
rest wavelength of the blend is 1549.479 A. 
g Blended with O ill] A1666 emission. 
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ally broad, with FWHM ~ 850 km s^ 1 , compared to 
the average of ~ 650 km s _1 seen in z ~ 2-3 galax- 
ies (Ste idel et al]l2010h . The reason for the unusually 
large line width is clear from the line profile shown in 
the right panel of Figure [SJ while typical Lya profiles of 
high redshift galaxies are purely redshifted and show a 
sharp blue edge, the Lya emission in BX418 shows a sig- 
nificant blueshifted component as well. The total range 
of velocities spanned by the line is large, ranging from ap- 
proximately —800 to +1500 km s" 1 . The line centroid 
is redshifted with Av = 300 km s _1 , somewhat lower 
than the z ~ 2 average of ~ 450 km s _1 ; this lower than 
average value is clearly due to the significant blueshifted 
emission, rather than to a lower characteristic outflow 
velocity. 

Given the ubiquitous nature of outflows in high red- 
shift galaxies, a great deal of energy has been devoted to 
modeling the radiative transfer of Lya photons in an ex- 
panding medium. This is a complex problem, and even 
the simplest models suffer from significant degeneracies 
in recovering a unique kinematic model from a given line 
profile. Given these degeneracies, and the fact that our 
resolution is too low to resolve any detailed structure 
that may be present in the profile, we do not attempt 
rigorous modeling. Nevertheless it is interesting to con- 
sider BX418 i n the context of these models . According 
to the work of iVerhamme et all (j2006l 12008ft . if the out- 
flow is considered to be a shell with velocity V cxp , the 
redshifted peak of the Lya emission falls at v ~ 2 x V exp ; 
this represents emission backscattered off the far side of 
the receding shell. There may also be secondary emission 
peaks at v ~ and v < — V CKp , corresponding to photons 
escaping from the blue and red wings of the approach- 
ing shell. In a more realistic model with a broad range of 
outflow velocities, the emission will span a broader range, 
with peaks at the velocities of lowest optical depth. The 
significant blueshifted emission in BX418 therefore sug- 
gests that the optical depth to Lya photons is unusually 
low in the foreground outflowing material. This is con- 
sistent with our conclusion in Section 17.1.11 above that 
the covering fraction of neutral gas is low, probably be- 
cause a higher than average fraction of the outflow is 
highly ionized. The low dust content undoubtedly also 
contributes to the apparent ease with which Lya photons 
e scape the outf l owing gas. 

ISteidel et al.l (|2010l ) have recently applied simple ana- 
lytic models to the kinematics of outflowing gas, and are 
able to reproduce observed interstellar absorption and 
Lya emission profiles with two gas components, one out- 
flowing and one at or near the systemic velocity of the 
galaxy. These results suggest that the observed line pro- 
files may be more dependent on bulk gas velocities than 
on detailed radiative transfer effects. We have not at- 
tempted to use these schematic models to reproduce the 
Lya line profile of BX418 in detail, but the models do 
demonstrate that significant blueshifted emission can be 
produced if the optical depth in the outflow is sufficiently 
low. 

Although BX418 is a sample of onl y one, it is inter- 
esting to compare with the results of IVerhamme et all 
(2008), who model the Lya emission in a sample of 11 
galaxies at 2.8 < z < 5. They suggest that the equiv- 
alent width and FWHM of Lya emission are correlated 
with the column density of neutral hydrogen iVm in op- 



posite senses, and that therefore systems with both high 
equivalent width and large FWHM should not be ob- 
served. They note that observed Lya emitters (LAEs, 
with W Lya > 20 A) have FWHM less than - 500 
km s _1 . BX418 clearly does not follow this trend. With 
WLya = 54 A, it is well above the threshold to be selected 
as an LAE, and yet it has F WHM ~ 850 km s" 1 . The 
IVerhamme et al.l ([2006, 2008) expanding shell models are 
unable to fit broad lines with the low column density 
required to produce a high equivalent width; however, 
these models assume unity covering fraction, which is 
unlikely to be the case in BX418, and it is probably the 
patchy distribution of optically thick gas which allows 
significant escape of Lya photons from the blueshifted 
approaching outflow, thus adding to both the equiva- 
lent width and FWHM of the line. This situation is 
unusual, but BX418 is not the only example; as noted in 
Section 18.31 galaxies with significant Lyman continuum 
emission also show blueshifted Lya emission (Steidel et 
al., in preparation). Models with a more realistic treat- 
ment of gas covering fraction may be more successful in 
reproducing the observed diversity of Lya profiles. 

Given this strong Lya emission and BX418's low dust 
content, we might expect that a large fraction of the Lya 
photons are able to escape the galaxy. We assess this by 
comparing the Lya and Ha fluxe s: for case B recombi- 
natio n, we expect Lya/Ha = 8.3 (|Ferland fe Osterbrockl 
|1985|) . We observe Lya/Ha = 3.7, a factor of ~ 2 lower 
than the case B predict ion, or escape fraction f e = 0.4. 
IVerhamme et a l. (2008) suggest that dust extinction is 
the primary factor in determining / e , finding that their 
sample is well fit by 

fe = 10 -7.71XE(B-V)_ (2) 

The correlation between dust co ntent and escape fr action 
is also seen in local galaxies by lAtek et all (l2009f). an d 
in a larger high redshift sample bv iKor nei et all (120101V 
who fi nd that the data scatter about the lVerhamme et afl 
( 2008) model, with the majority of points having a lower 
escape fraction tha n predicted by Equa tion [2] (all calcu- 
lations assume the iCalzetti et al.l 120001 extinction law) . 
From this relation and our inferred E(B — V) = 0.02 we 
would predict /«, = 0.7, higher than the observed value; 
like most of thelkornei et al.) (120101) sample, BX418 falls 
below the IVerhamme et alJ()2008D prediction. However, 
given the significant uncertainties in E{B — V), the pre- 
dicted f e is within our uncertainties. This comparison 
also assumes that Ha and Lya emission have the same 
spatial extent, which may not be the case because Lya is 
subject to greater scattering. Our observed f e is there- 
fore a lower limit, since the spectroscopic slit may be 
missing a larger fraction of the Lya photons. 

7.1.3. Intervening Absorption and the Physical Extent of 

Outflows 

One of the stronger absorption features in the spec- 
trum, with equivalent width W = 2 A, falls at 4926 A 
(see Figure [7]) . There is no intrinsic absorption feature 
at this wavelength in the rest frame of BX418, as shown 
by the composite spectrumQ We identify this line, and 

12 This absorption feature prev ents the detection of any possible 
N iv] A1487 emission; see Section EH 
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a corresponding feature at 3859 A, as intervening C IV 
and Lya absorption from the foreground galaxy Q2343- 
BX442, with z = 2 176. BX442, identified in the Ha 
survey of lErb et al.l (|2006d) . is separated from BX418 by 
an impact parameter of 17", or 140 kpc. Other galaxies 
we find near BX418 are known or likely to be at lower 
redshifts, given their colors and magnitudes. 

Absorption systems along the line of sight to back- 
ground galaxies and QSOs can be u sed to assess the 
spati al extent of galactic o utflow s (e.g.. lAdelberger et al.l 
2003t lAdelberger et all [2005t ISimcoe et all 120061: 
Steidel et al.l I201QT K Close pairs of background QSOs 
and foreground galaxies are rare but offer high S/N; 
galaxy-galaxy pairs are much more common but are 
more diffi cult because o f the fa intness of the background 
object. Steidel et al.l (|2010D have recently stacked 
spectra of galaxy-galaxy pairs as a function of impact 
parameter, to determine the average absorption of 
various transitions with distance; they find significant 
absorption in both C IV and Lya at an impact pa- 
rameter of ~ 60 kpc, and at ~ 100 kpc the average 
spectrum still shows clear absorption in Lya. Individual 
pair s can show much stron ger a bsorption; for exam- 
ple, lAdelberger et al.l (|2005[ ) and ISimcoe et all (|2006D 
identify strong absorption from C IV, Lya, O VI and 
other ions in a QSO spectrum, and associate it with a 
massive galaxy at an impact parameter of 115 kpc. With 
equivalent width ~ 2 A, the absorption from BX442 
is similarly well above the average value (absorption 
of this strength would hav e been easily detec table in 
the composite spectra of Steidel e t al.l l2010f ). The 
sample of individual pairs of galaxies with detected 
foreground absorption is still too small to test for any 
correlations between the properties of the foreground 
galaxies and the strength of the observed absorption, 
but we note that BX442 is among the oldest and most 
massive galaxies in the sample, with an estimated age 
of 2.8 Gyr and stellar mass M* = 1.1 x 10 11 M ; it 
has therefore had considerable time in which to drive 
material to large distances. Given BX442's large stellar 
mass, it probably occupies a dark matter halo more 
massive than the z ~ 2 average, (Mhaio) = 9 x 10 11 M Q 
(jConrov et al.ll2008D ; for a halo of mass - 2 x 10 12 M , 
the virial radius is r vlr ~ 125 kpc (see calculations by 
iSteidel et all 120101 ) . Thus it appears that BX442 has 
driven significant quantities of gas to the virial radius 
and beyond. 

7.2. He II Emission 

One of the most striking differences between the UV 
spectrum of BX418 and the average z ~ 2 galaxy spec- 
trum is BX418's strong He II A1640 emission, which has 
equivalent width Wq = 2.7 A and flux relative to H/3 
-FHeii/-FH0 ~ 0.3. The line is shown in Figure [5J with 
the z ~ 2 composite spectrum overplottcd for compar- 
ison. Broad He II A1640 emission is characteristic of 
high redshift galaxies; the weak emission in the compos- 
ite shown in Figure [9] has equivalent width Wq ~ 0.6 
A, and He II emission is clearly seen in the composite 
spectrum of z ~ 3 Lyman break galaxies analyzed by 
IShaplev et all (f2003h. with Wq = 1.3 ± 0.3 A (as mea- 
sured by IBrinchmann et~aHl2008bQ and FWHM - 1500 
km s _1 . Because of its broad width, the line is at- 



tributed to the fast, dense stellar winds of Wolf-Rayet 
stars. Broad He II emission with Wq = 2.45 ± 0.22 A is 
also seen in the spectrum of the z = 2.7 lensed galaxy 
know n as the 8 o'clock arc (|Dessauges- Zavadsk v et all 
I2010D , and at similar strength i n the spectrum of th e 
lensed z = 3.8 galaxy observed bv lCabanac et al.l (|2008ft . 
an object which otherwise appears similar to cB58. 

While the He II emission in BX418 is clearly broad, it 
also has a narrow peak not seen in the composite spectra 
or in the high resolution spectrum of the 8 o'clock arc. 
As shown in the right panel of Figure [SJ the line is well 
fit by the superposition of two Gaussian components: a 
broad line with FWHM ~ 1000 km s _1 accounts for 
~ 75% of the line flux, while a narrow, unresolved com- 
ponent makes up the remaining ~ 25%. Compared to a 
single Gaussian, the superposition of these two compo- 
nents provides a better fit to the line at the 99% con- 
fidence level; a single Gaussian underpredicts both the 
broad wings and the narrow peak of the line. As with 
the composite spectra, the width of the broad component 
indicates an origin in the winds of W-R stars, while the 
narrow component is likely to be nebular emission. We 
discuss the stellar and nebular He II emission separately 
below. 

7.2.1. Stellar He II Emission 

The strength of the stellar He II line undoubtedly de- 
pends on the number of W-R stars producing the emis- 
sion, and probably depends on the type and metallicity of 
these stars as well. Because the W-R phase is very short, 
age and star formation history are also crucial. There 
has been considerable recent work addressing these is- 
sues, both locally and at high redshift, but progress is 
hampered by the lack of local observations of individual 
W-R stars of low metallicity. 

Although detailed observations are still needed, cur- 
rent models indicate that line luminosities of W-R 
stars decrease with decreasing met allicity, due t o the 
reduced density of stellar winds (S mith et al.l 120021 : 
lCrowther~fe Hadficld 2006). Studies of lensed galax- 
ies and composite spectra are broadly consistent with 
this trend, as the only strong stellar He II emission 
seen in a lensed galaxy with a determination of metal- 
licity belongs to the massive and m etal-rich 8 o'clock 
arc ()Dessauges-Zavadskv et al.l 120101 : metallicity mea- 
surements are not a vailable for the galaxy observed by 
iCabanac et aLll2008h . This emission is also stronger than 
than seen in the z ~ 3 composite spectrum, which prob- 
ably has a lower average metallicity. 

The pre cise form of the line scaling with metallicity is 
unknown. Brinch mann et all (|200 8b) combine models of 
He II A1640 luminosities with a variety of star formation 
historie s in an effort t o acc ount for the He II emission 
seen by IShaplev et al.l (|2003| ): they find that, while He II 
emission peaks a few Myr after a burst of star forma- 
tion, even for more extended star formation histories the 
line should be detectable with an equivalent width of 
1-2 A, at metallicities above Z ~ 0.5 Z©. While this 
model may account for the He II emission seen in the 
z *~ 3 and z ~ 2 composite spectra and in the 8 o'clock 
arc, it is more problematic in the case of BX418, with 
Z = 0.2 Zq. At this low metallicity, the He II emission is 
never predicted to exceed Wq ~ 0.5 A, even immediately 
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Figure 9. Left: The strong He II A1640 emission line (black), compared to the broad and far weaker emission in the composite spectrum 
of 966 z ~ 2 galaxies (red). Right: The line is well fit by a superposition of two Gaussian components (solid red line). The broad component 
(long-dashed dark blue line) has FWHM ~ 1000 km s — 1 , while the narrow component (short-dashed green line) is unresolved. We attribute 
the broad emission to W-R stellar winds, and the narrow component to nebular He II emission. 



after a burst of star formation. More recent stellar pop- 
ulation synthesis models w hich include binary evolution 
(|Eldridge fc Stanwavll2009D offer better agreement with 
our observations, however. These models so far include 
only instantaneous bursts rather than more extended star 
formation histories, but they find that with binary evolu- 
tion taken into account, W-R features are both stronger 
and occur at later times. For a metallicity Z = 0.004, 
close to the metallicity we find for BX418, their models 
predict an equivalent width Wq — 3 A at an age of ~ 10 
Myr. The equivalent width is stronger at lower metal- 
licities because the W-R phase lasts longer; WNL stars 
form late, through binary evolution, and have a longer 
lifetime at lower metallicities because it takes longer for 
the weak stellar winds to remove the stars' hydrogen en- 
velopes. Thus the lower metallicity stars make a greater 
contribution to the integrated spectrum. 

While these models may account for the He II emission 
seen in BX418, the situation remains complicated. Com- 
parison of BX418 with the z ~ 2 composite spectrum 
supports the scenario in which low metallicity stars make 
a greater contribution to the He II equivalent width. We 
infer from the z ~ 2 mass-metallicity relation that the 
average metallicity of the ga laxies making up the z ~ 2 
composite is about .5 (|Erb et al.l 12006a!) . and the 
lEldridge fc Stanwavl (|2009f) models predict significantly 
weaker He II emission at this metallicity. The observa- 
tions of lensed galaxies described above suggest the op- 
posite picture, however, since the only detection of He II 
emission is in the most metal-rich object observed. In ad- 
dition, the effect of varying star formation histories still 
needs to be taken into account; the instantaneous burst 
model is certainly not appropriate for typical galaxies in 
the z ~ 2 sample and is probably not the best model for 
BX418 either. Even for BX418, it remains to be de- 
termined how much of the He II line strength is due 
to metallicity, and how much may be due to observa- 
tion at a particular, brief stage of evolution. Both addi- 
tional spectra of low metallicity galaxies and models of 



the line strength for a variety of star formation histories 
are needed to address this question. 

7.2.2. Nebular He II Emission 

Decomposition of the He II line profile indicates that 
about 25% of the flux arises from unresolved nebu- 
lar emission, for which there is no evidence in the 
composite spectra at either z ~ 2 or z ~ 3. This 
is an immediate indication of an unusually hard ion- 
izing spectrum, with significant numbers of photons 
with energies greater than 54.4 eV. Nebular He II 
emission is occasionally observed in local W-R galax- 
ies, and app ears to occur primari l y in lower metallic- 
ity systems ([Campbell et al.l 119861: iGarnett et al.l 119911 : 
iChandar et alj|2004t iBrinchmann et al.ll2008aD . Even lo- 
cally, however, explanations for the source of the ionizing 
photons differ . Models of the ionizing fluxes of W-R stars 
([Smith, et aTl2"0 02) show that flux above 54 eV is a strong 
function of metallicity; in the weaker, less dense winds 
of metal-poor W-R stars, He ++ does not recombine as it 
does in a denser wind, and significant flux i s emitted at 
energi es above the He + edge. The models of lSmith et al.l 
(2002) also predict that nebular He II emission will be 
seen only immediately after an instantaneous burst of 
star formation; for continuous star formation, the frac- 
tion of sufficiently energetic photons is so small that neb- 
ular He II emiss ion is not expected to be detectable at 
any metallicity. IBrinchmann et all ([2008af ) also investi- 
gate the origin of nebular He II emission, using a large 
sample of galaxies showing W-R features from the Sloan 
Digital Sky Survey. Using timescales derived from H/3 
equivalent widths, they argue that nebular He II emis- 
sion arises earlier than would be expected if its primary 
source were W-R stars, and that therefore low metallic- 
ity O stars must produce significant numbers of photons 
above 54 eV. At least qualitatively speaking, either of 
these explanations can probably account for the nebular 
emission in BX418. Its low metallicity should allow the 
production of He-ionizing photons from low metallicity 
stars, and the unusual strength of the stellar He II emis- 
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Figure 10. Predicted nebular He II emission as a function of metallicity and ionization parameter, from Cloudy photoionization modeling. 
The colors represent different metallicities as labeled at upper left, and the shaded region shows the estimated He Il/H/3 ratio of BX418, 
including uncertainties. For a metallicity Z ~ 0.2 Zq, an ionization parameter log U ~ — 1 is required to match observations. 



sion suggests that we are observing a strong burst of star 
formation. 

We further investigate the conditions needed to pro- 
duce significant nebular He II emission with photoion- 
ization models from Cloudy. As described in Section [SJ 
we use an input ionizing spectru m from Starburst99 , in- 
corporating the ionizing fluxes of lSmith et "all (|2002ft and 
using a metallicity Z = 0.004. We caution that in inter- 
preting the Cloudy results, it is important to remember 
that the model output is only as reliable as the input, 
and the input ionizing spectru m above 54.4 eV s hould 
be considered highly uncertain. iSmith et al.l (|2002f ) note 
that the flux above 54.4 eV depends sensitively on the 
wind density and on the form of the law used to scale 
the winds with metallicity; this scaling law isn't known 
precisely even for O stars, and for W-R stars it is not 
known at all. They also point out that the grid of mod- 
els excludes the rare, hot W-R stars with weak winds 
that have significant ionizing fluxes above 54 eV; these 
are very rare in the Galaxy and LMC, but could be more 
significant at lower metallicities. It is clear that a more 



thorough understanding of the source of these energetic 
photons locally is required in order to interpret results 
at high redshift. 

With these cautions in mind, we proceed with the mod- 
eling assuming that 25% of the He II emission is nebu- 
lar, for a flux relative to H/3 of log (-Fh g ii/-Fh^) — —1-1. 
The Cloudy models then predict the nebular He II flux 
as a function of metallicity and ionization parameter. 
The results of the modeling are shown in Figure [TUJ 
the He II/H/3 ratio increases with both metallicity and 
ionization parameter, and for the metallicity of BX418, 
an extremely high ionization parameter logU ~ — 1 is 
required to account for the observed emission. Inter- 
estingly, the increase of He II/H/3 with metallicity and 
ionization parameter is not primarily driven by an in- 
crease in He II luminosity. Rather, at very high ion- 
ization parameters (log U > —2) and in the presence 
of dust grains, hydrogen-ionizing photons are selectively 
depleted by energy-specific grain opacities peaking near 
912 A, resulting in a decreasing H ^ flux with both ion - 
ization parameter and metallicity (|Bottorff et alj fl998). 
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The metallicity dependence arises because higher metal- 
licity models are assumed to also have more grains. The 
implications of this effect are discussed further in Sec- 
tion 18. 1| for now, we note that entirely dust-free mod- 
els do not successfully reproduce the observed line ratios 
of BX418, and that in spite of the inferred low extinc- 
tion such models seem unrealistic given the metallicity of 
~ 1/6 7iq and the known correlation between metallicity 
and dust content. Also recall that log U > — 2 is required 
by the lower limit on the [O III]/[0 II] ratio (Section [GT]). 

Log U ~ — 1 is much higher than the ionization pa- 
rameters inferred for more typical z ~ 2 galaxies (which 
are in turn higher than thos e inferred for l ocal star- 
forming galaxies); for example. lHainline et al.l ()2009l ) use 
the [O III] / [O II] ratio of two lensed and apparently rep- 
resentative z ~ 2 galaxies to infer \ogU ~ —2.5. While 
unusual, this high ionization parameter is not inconsis- 
tent with other observations, and may be expected from a 
compact and energetic ionizing source such as a cluster of 
massive, low metallicity stars. We discuss this high ion- 
ization parameter further in Section I8.1[ and emphasize 
in conclusion that regardless of the uncertainties associ- 
ated with the input ionizing spectrum and the photoion- 
ization modeling, the fact remains that the He II emission 
seen in BX418 is exceptionally strong and indicative of 
an ionizing spectrum that must be significantly harder 
than that seen in more typical z ~ 2-3 galaxies. 

7.3. Nebular Emission Lines, Direct Metallicity and 
C/O Abundance 

The UV spectrum of BX418 is also notable for its 
strong nebular emission from O III] AA1661, 1666 A and 
C III] AA1907, 1909 A. Inspection of Figure shows that 
the C III] emission is far stronger than that observed 
in the composite z ~ 2 spectrum, while O III] emission 
is not present in the composite at all (it is, however , 
detected in the z ~ 3 composite of iShaplev et al.lfeOOSL 
and is noticeably stronger in their subset of galaxies with 
strong Lya emission; see Section l5T2")) . We use the Cloudy 
models to investigate the conditions under which such 
strong O III] and C III] emission arises. 

The photoionization models show that O III] emission 
increases strongly with ionization parameter, and de- 
creases with metallicity; at Z = 0.5 Z Q and log U = —2.5, 
the combined strength of the O III] doublet is roughly 
~ 5% of the flux of H/3, while for Z = 0.2 Z Q it reaches 
the observed flux of ~ 25-30% of H/3 for ionization pa- 
rameters ranging from log U ~ —1.6 to —0.8. This is only 
partially dependent on the decrease in H/3 at very high 
ionization parameters discussed above; much of the in- 
crease in the O III] flux with ionization parameter occurs 
well before the photon depletion effect becomes signifi- 
cant. This result is consistent with the high ionization 
parameter inferred from the He II emission, and accounts 
for the lack of O III] emission in the z ~ 2 composite as 
well. We also note that the detection of the O III] line is 
complicated by the presence of Al II absorption at 1670 
A; the line is undoubtedly more easily detected when the 
low ionization absorption lines are weak, as in BX418, or 
in higher resolution spectra. 

The ratio O III] AA1661, 1666/ [O III] A5007 is sensitive 
to the electron temperatu re. Using the calibration of 
iVillar-Martm et "ahl (|2004l . see their Figure 1), we find 



T e ~ 15, OOOIyo" K, in agreement with the best fitting 
photoionization models, which also have T e ~ 15, 000 K. 
This determination of the electron temperature allows a 
direct me asurement of the m etallicity via the T e method; 
following ILTotov et al.l (|2006f ). we find 12 + log(0/H) = 
7.8 ±0.1, in agreement with the i?23-derived metallicity 
of 12 + log(0/H) = 7.9 ± 0.2 (Section GO}. 

The modeling also shows that C III] emission is 
strongly dependent on both metallicity and ionization 
parameter, with the line strength increasing with ioniza- 
tion parameter at all metallicities, peaking at Z ~ 0.2 
Zq, and decreasing at both higher and lower metallici- 
ties. Thus BX418 appears to have the optimal conditions 
for strong C III] emission, with an extremely high ion- 
ization parameter and a metallicity that is low but not 
too low. Further investigation of the models shows that 
the metallicity dependence of C III] emission is largely a 
dependence on the abundance of carbon, which we study 
in more detail below. 

Other high ionization nebular emission lines are occa- 
sionally seen in high redshift gal axy spectra. Weak N IV] 
A1487 emission is seen in cB58 (jPettini et al.ll2000h . but 
because of the intervening absorption feature discussed 
in Section l"7. 1.31 we are unable to tell whether or not it is 
present in BX418. The Cloudy photoionization models 
predict weak N IV] emission, at the level of a few per- 
cent of H/3, depending on the N/O ratio. We have also 
checked for the presence of N III] A1750 emission, finding 
an upper limit of F(N III] /H/3) < 0.1. Unfortunately this 
is not sensitive enough to constrain the nitrogen abun- 
dance; for all metallicities and all values of N/O consid- 
ered, the Cloudy models predict F(N III] /H/3) < 0.1. We 
discuss the possibility of nebular C IV emission in Section 
1731 

7.3.1. The C/O Abundance Ratio and the Origin of Carbon 

The origin of carbon is a subject of some contro- 
versy. It is produced primarily through He burning via 
the triple a reaction, which may occur in both massive 
(M > 8 M Q ) and low to intermediate mass (1 < M < 8 
M©) stars; however, there is considerable disagreement 
over the mass range of the stars responsible for most 
carbon production. The C/O ratio is observed to in- 
crease with increasing O/H in both individual stars and 
in galaxies and H II regions. This trend has been stud- 
ied by many authors, using a variety of chemical evolu- 
tion models and stellar yields. Some find that the data 
are best explained by C arising almost exclusively from 
massive stars, in which case the trend in C/O vs. O/H 
is due to metallicity-dependent stellar winds, as mass 
loss and ISM enrichmen t are gre ater at higher metallic- 
ities (jHenrv et al.l [20001. Others (jChiappini et al.l 120031 ) 
reach the opposite conclusion, finding that carbon pro- 
duction is dominated by low to intermediate mass (1 < 
M < 3 M Q ) stars, and that the C/O vs. O/H trend is 
therefore an evolutionary effect due to the delayed re- 
lease of carbon relative to oxygen (which is produced 
almost exclusively by m assive stars) in yo unger and 
less metal-rich systems. iCarigi et al.l (|2005| ) find that 
massive and low to intermediate mass stars contribute 
rough ly equal amoun t s of c arbon in the solar vicinity, 
while Akcrman et al. (2004) find that the metallicity- 
dependent mass loss of high mass stars is the primary fac- 
tor in producing the C/O vs. O/H trend, with some sec- 
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Figure 11. The ionization states of oxygen (left panel) and carbon (middle panel) as a function of ionization parameter, from the Cloudy 
models. The right panel shows the ionization correction factor (see EquationUJl as a fu nctio n of ionization parameter. The colors represent 
mctallicities ranging from Z = 0.05 Zq (blue) to Z = 0.5 Z0(red), as labeled in Figure [Tol 



ondary contribution from delayed C production in lower 
mass stars. 

Measurement of the C/O ratio in a very young star- 
burst can provide valuable discrimination between these 
models, since any time delay in carbon production will 
be most pronounced at early times and with high star 
formation rates. We estimate the C/O abundance of 
BX418 via the ratio of O III] AA1661, 1666 A to 
C HI] AA1907, 1909 A follo wing the method des cribed 
bv lGarnett et all £[995) and IShaplev et all ([20031 ) . who 
derive the C +2 /0 +2 ionic abundance ratio from the 
emission-rate coefficients for collisionally excited emis- 
sion lines: 



C+ 2 

0+2 



= 0.15 e 



-1.106V* J (° E] AA1907 > 1909 ) 
1(0 in] AA1661,1666)' 



(3) 



where t — T e /10 4 K. Note that t his ex pression differs 
slightly from that of Garnctt et all (fl9 95'). which includes 
only O III] A1666; we include both members of the O III] 
doublet. This ionic abundance ratio will be equal to the 
true C/O abundance only if C and O have the same 
relative ionization. This may be approximately true in 
many circumstances, but is not necessarily the case. + 
and +2 have higher ionization potentials than C + and 
C+ 2 (35.1 and 54.9 eV vs. 24.4 and 47.9 eV respectively), 
so in regions ionized by a hard ionizing spectrum, as is 
apparently the case in BX418, a significant amount of 
carbon may be in the form of C +3 , and the C +2 /0 +2 
ionic abundance ratio would then underestimate the true 
C/O abundance. In order to correct for this effect it is 
necessary to ap ply an ionization corr ection factor (ICF) 
as described bv lGarnett et all ([19950 : 



C 
O 



c+ 2 

0+2 
0+2 
0+2 



X(C 



+2^1 



Y(0+ 2 ) 
x ICF, 



(4) 



where A(C+ 2 ) and X(0 +2 ) are the C+ 2 and 0+ 2 volume 
fractions, respectively. 

We use the Cloudy models to examine the ionization 
states of C and O, and estimate the ionization correc- 
tion factor as a function of ionization parameter. The 
results are shown in Figure [TTJ The two left panels show 
the ionization fractions of O and C as a function of ion- 
ization parameter; the difference due to oxygen's higher 



ionization potential is apparent, as for log U > —2, nearly 
all the oxygen is in the form of +2 , while the fraction 
of C +2 declines at \ogU > —2, becoming equal to the 
fraction of C +3 around logC/ ~ — 1. We translate these 
fractions into the ionization correction factor in the right 
panel of Figure 111! For moderate ionization parameters 
around \ogU ~ — 2, the ICF ~ 1 (with little dependence 
on metallicity) , and C +2 /0 +2 will be a good approxi- 
mation for C/O. Indeed, i n the sample of dwarf galaxies 
studied bv lGarnett et all JT995), the range of ICF found 
is only 1.06-1.33. On the other hand, with increasing 
ionization parameter the ICF becomes more significant, 
and the metallicity dependence increases as well. 

There is considerable uncertainty associated with the 
derivation of the ICF for BX418. Our primary constraint 
on the ionization parameter comes from the He II flux; 
this determination of log U relies heavily on the un- 
certain input ionizing spectrum used by the models, as 
discussed in Section 17.21 For Z — 0.2 Z Q , the best- 
fitting models require approximately log U = — 1.0±0.15. 
The corresponding ionization correction factor is there- 
fore ICF = 1.7 ± 0.2. For the electron temperature we 
adopt T e = 15,0001^ K as indicated by the O III] 
AA1661, 1666/[0 III] A5007 ratio. Including these uncer- 
tainties in the ICF, in T e , and the flux uncertainties in 
the C III] and O III] lines, we find C/O = 0.24 ± 0.05, 



and log (C/O) = -0.62 



+0.08 
-0.1 



This is in agreement with 



the input C/O ratio of the best-fitting Cloudy models, 
all of which have C/O = 0.19. 

We show the log (C/O) vs. log(0/H) diagram in Fig- 
ure 1121 using the oxygen abundance determined by the 
electron temperature method for consistency with local 
observations. BX418 is shown by the large magenta cir- 
cle, and the remaining symbols show local galaxies as de- 
tailed in the caption. The shaded region shows the range 
of C/O determined fro m the LBG composite spectrum 
bv lShaplev et~aT1 (I2003T) . log (C/O) = -0.68±0.13; there 
is no reliable O/H measurement for this set of galaxies, 
and n o ionization correction was applied. IShaplev et all 
(2003) also found a somewhat lower value of C/O for 
their subset of galaxies of galaxies with strong Lya 
emission, log (C/O) = —0.74 ± 0.14, again with no 
ionization correction applied (for comparison, we find 
log(C +2 /0 +2 ) = —0.87 without the ionization correc- 
tion for BX418). 

BX418 shows excellent agreement with the local trend 
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Figure 12. The carbon-to-oxygen ratio vs. the oxygen abun- 
dance. BX418 is shown by the large magenta circle. The other 
points show values for local g alaxies. Blue triangles ar e lo- 
cal dwarf irregular gala xies from Garnctt et al. (1995, 1993) an d 
IKobuln ickv & Skillmarj] HMB ), a nd the open light and dark pur- 
ple circles show local spirals from Garnctt ct al. (1999), with two 
different extinction corrections. The black star shows the default 
value of the H II region abundance set used by Cloudy, correspond- 
ing to the Orion nebula; the most recent estimate of this ratio in 
Orion is log (C/O) = -0.21 ( EstebaneFaD EH) . The shaded 
region shows the value determined from the z ~ 3 LBG composite 
spectrum by Shaplcv ct al. (2003). 

in C/O vs. O/H, in spite of its differences with respect 
to the local sample. The local objects it most resem- 
bles in abundance ratios are H II region s of dwarf emis- 
sion line galaxies (G arnett et al.l [19951) . While a de- 
tailed comparison with such objects is difficult, local 
H II galaxies probably have more extended star forma- 
tion histories than BX4I8, extending to at least 1 Gyr 
and consisting of short bursts alt ernating with quies- 
cent periods frerlevich et all 120041: iWestera et all 120041: 
IPerez-Montero et al.1 ^OloT we assess the possibility of 
an earlier episode of star formation in BX418 in Sec- 
tion 14.11 and find that it is a poorer fit to the observed 
SED than a single stellar population, but we cannot rule 
out this possibility entirely). The star formation activ- 
ity of these galaxies is also less extreme than that seen 
in BX 418, as indicated by their lower ionization param- 
eters (iPerez- Montero fc Diaz 2003; Pcrcz -Montero et al.l 
I2010T) . If a delay in carbon production by low mass stars 
were the primary cause of the increase in C/O with O/H, 
we would likely see a significantly lower C/O ratio in a 
very young, extreme starburst. Furthermore, the best- 
fit age for BX4I8, ~ 38 Myr, allows only massive stars 
with M > 10 Mq to have left the main sequence; an 
age of 100 Myr allows for carbon enrichment from stars 
with M > 6 M . It is therefore unlikely that carbon 
production i s dominated by 1 < M < 3 M Q stars as 
sugg ested by iChiappini et all (|2003l ). Rather, we concur 
with lAkerman et al.l ( 20041) . who conclude that the ap- 
parent universality of the metallicity dependence of C/O 
in objects of widely varying star formation histories is 
more likely to be due to variable stellar yields than to 
evolutionary effects. 

7.4. Stellar Wind Lines 



The most prominent stellar features in UV spectra of 
high redshift galaxies are generally the wind lines from 
massive stars; the high ionization lines of N V A1240, 
Si IV AA1394, 1402 and C IV AA1548, 1551 show distinc- 
tive P Cygni pr ofiles indicative of sta rs with masses of 
at least 30 M (jLeitherer et alJ[l995| ). These lines de- 
pend sensitively on the properties of the massive stellar 
population, so in principle they may provide a great deal 
of information. In practice, however, their interpretation 
is difficult. The strength of these features depends on 
metallicity, the age of the starburst, the slope and upper 
mass cutoff of the IMF, and possibly on time -dependent 
dust extinction as well (Lcithcrcr ct al. 2002). The lines 
themselves are complex blends of stellar wind and pho- 
tospheric features and interstellar absorption; nebular 
emission may also contribute to the line profiles. The 
situation is further hampered by our lack of detailed un- 
derstanding of the scaling of stellar winds with metallic- 
ity, and the limited number of observations of massive 
stars at low metallicities. Given these substantial uncer- 
tainties, and the relatively low resolution of the spectrum 
which precludes accurate separation of the stellar and in- 
terstellar lines, we cannot draw robust conclusions from 
the stellar wind features in the UV spectrum of BX418. 
Nevertheless these lines show intriguing differences from 
typical high redshift galaxy spectra, as we show below. 

Figure[felshows the N V A1240, Si IV AA1394, 1402 and 
C IV AA1548, 1551 lines in the spectrum of BX418, with 
the composite spectrum of ~ 1000 z ~ 2 galaxies over- 
plotted (in red) for comparison. For all three lines, the 
most notable difference between BX418 and the compos- 
ite spectrum is the strength of P Cygni emission, which 
is stronger in BX418 in each case (the emission occurs at 
1244, 1405 and 1553 A in the N V, Si IV and C IV lines 
respectively). The broad, blueshifted absorption, most 
apparent in N V and C IV, is very similar in BX418 and 
the composite (the differences in the narrower interstel- 
lar abs orption lines of Si IV and C IV are discussed in 
Section l7XTj) . 

Broadly speaking, the behavior of these lines with vary- 
ing stellar population properties is understood: the P 
Cygni emission is strongest in a young stellar population 
of 3-5 Myr and stabilizes by ~ 20 Myr, increases with any 
variation of the IMF that increases the number of massive 
stars (a flatter slope or higher upper mass cutoff), and de- 
creases at lower metallicities. The strength of absorption 
depends primarily on the density of the stellar wind and 
therefore on the mass loss rate. Mass loss rates are lower 
at lower metallicities, and thus the absorption troughs 
are weaker at lower metallicities. The absorption profiles 
also show a mild age dependence, particularly in the case 
of C IV, with stronger absorption in younger starbursts. 
However, significant work remains to be done in quan- 
tifying these dependencies and in creating a wide range 
of models for compa rison with observatio n. The flexible 
Starburst99 models (jLeitherer et al.| [l999) allow the cre- 
ation of synthetic UV spectra with a variety of ages and 
IMFs, but include only two metallicities, solar and the 
averag e LMC/SMC me tallicity of Z ~ 0.4 Z . More re- 
cently, iRix et all ((2004) combined Starburst99 with the 
model atmosphere code WM-BASIC to create a suite of 
synthetic spectra covering a wider range of metallicities; 
however, this work was focused primarily on the stellar 
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Figure 13. Profiles of the stellar wind lines N V A1240 (left panel), Si IV AA1393, 1402 (middle panel), and C IV AA1548, 1551 (right 
panel). The spectrum of BX418 is shown in black, and the composite spectrum of ~ 1000 z ~ 2 galaxies is shown in red for comparison. 



photospheric lines, and the authors caution that a more 
detailed study of the wind lines is still needed. Because of 
the focus on photospheric lines, processes such as shock 
emission that affect only the high ionization wind lines 
were not included in the models. These limitations have 
been addressed to some extent in the latest g eneration of 
the WM-BASIC code (jLeitherer et all 120101 ) . However, 
the wind line profiles generated by running Starburst99 
with this latest version of WM-BASIC do not differ ap- 
preciably from those used in the present work, and none 
of the discussion below is altered by the use of these lat- 
est models. 

We compare the profiles of BX418's stellar wind lines 
with various models in Figure [Ml The top three rows are 
Starburst99 LMC/SMC metallicity models with Z ~ 0.4 
Z ; this is a factor of ~ 2 higher than the metallicity 
of BX418, and these models may be expected to have 
stronger emission and absorption than comparable mod- 
els of lower metallicity. Nevertheless they are illustrative 
of the behavior of the wind lines with variations in age 
and IMF. The top panels of Figure (dark red lines) 
show the standard model which has been shown previ- 
ously to give generally good agreement with high redshift 
galaxy spectra; this model has a Salpcter slope and an 
upper mass cutoff of 100 M , and an age of 50 Myr (other 
studies have assumed an age of 100 Myr; there is almost 
no difference between the two models). The second set 
of panels (blue lines) show a model with the same IMF 
but a younger age of 5 Myr, and the third row (green 
lines) shows a 50 Myr old spectrum with a flatter slope 
of a = 1.4, resulting in more massive stars. The primary 
effect of these two variations is to strongly increase the 
predicted emission. The bottom row shows an interpo- 
la tion between the Z = 0.05 Z and Z = 0.2 Z models 
of IRix et aLI 1)20041 ) . to create a Z ~ 0.1 Z© model in 
which the depth of the broad C IV absorption profile ap- 
proximately matches observations. Note that this model 
predicts no emission, except in the N V line. As we 
discuss further below, none of the models considered si- 
multaneously match the strengths of the absorption and 
emission components of the P Cygni features. 

As noted above and shown in Figure Q21 the depth of 
absorption is almost identical in BX418 and in the z ~ 2 
composite spectrum; this is surprising given that the lat- 
ter almost certainly has a significantly higher average 
metallicity. The absorption profiles do also depend on 
age, but models indicate that the strength of the absorp- 



tion is much more sensitive to metallicity than to age. 
It is possible that differences in the average age of the 
stellar populations cancel out the differences in metallic- 
ity to create nearly identical absorption profiles, but this 
would be a surprising coincidenc e . Mor e likely, the fact 
that the Z ~ 0.1 Z IRix et all (|2004[ ) model matches 
the composite spectrum as well as BX418 is an indica- 
tion that our understanding of the age and metallicity 
dependence of these lines remains incomplete. 

The Z ~ 0.1 Zq model matching the depth of the P 
Cygni absorption predicts no emission in either the C IV 
or Si IV lines, clearly at odds with observations. Stan- 
dard models with a Salpeter slope and an age of at least 
50 Myr also underpredict the emission component of the 
lines. Models matching the strong emission require ei- 
ther a very young age (Figure 1141 blue model in second 
row), or an IMF with more massive stars (green model, 
third row; note that both of these models overestimate 
the metallicity of BX418, so even more extreme varia- 
tions may be required to match the observed emission). 
As discussed in Section ITT21 a young burst superimposed 
on a somewhat older population may be required to pro- 
duce the strong stellar He II emission; therefore a young 
age such as that seen in the second row of Figure [TJ] may 
not be unreasonable. Perhaps the simplest explanation, 
however, and the one most likely to account for both the 
absorption and emission, is that the C IV and Si IV emis- 
sion is at least in part nebular. The best fitting Cloudy 
models predict C IV emissio n at the level of ~ 8% of H/3 
and - 15% of C III] (see also lLeitherer et al.|[2002l ). This 
is consistent with observations, although we may expect 
nebular emission to be difficult to detect when superim- 
posed over the strong interstellar and stellar absorption; 
on the other hand, as resonance lines, these transitions 
are subject to multiple scatterings, and like Lya, they 
may be most likely to escape the galaxy when they are 
redshifted away from their resonant wavelength. This 
would effectively mimic the stellar P Cygni emission pro- 
file. No nebular N V em ission is predicted by the models, 
but the IRix et al.l (|2004D models predict stellar N V emis- 
sion even at metallicities when there is none from C IV 
or Si IV. 

Finally, we note that the failure of models to reproduce 
all of the features of the wind line profiles is not unique to 
BX418. Most previous studies of UV spectra of high red- 
shift galaxie s have noted either an overpredict i on of the 
absorption (jPettini et al.1 120001: IShaplev et ail I2003D or 
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Figure 14. A comparison of the profiles of the wind lines N V A1240 (left panels), Si IV AA1393, 1402 (middle panels), and C IV AA1548, 
1551 (right panels) with synthetic spectra. The top three rows are continuous star formation Starburst99 models using the composite 
LMC/SMC metallicity of Z ~ 0.4 Zq. The top row (dark red model) shows a 50 Myr old starburst and an IMF with Salpeter slope and 
upper mass cutoff of 100 Mq . The second row (blue) shows a 5 Myr old starburst and an IMF with Salpeter slope and upper mass cutoff 
of 100 Mq, a nd the third row (green) shows a 50 Myr old starburst with a flatter IMF slope of a = 1.4. The bottom row (orange) shows 
a model from Rix ct al. (2004) with Z ~ 0.1 Zq, an age of 50 Myr, and a Salpeter IMF; this model was chosen to match the depth of the 
broad absorption in the C IV profile. The spectrum of BX418 is shown in black in all panels. 



an underprediction of the emission (jQuider et alJ feOOQ). 
Metallicity effects and nebular emission have been in- 
voked to account for the discrepancies. Future observa- 
tions of massive stars over a wide range of metallicities 
and improved models of stellar winds across this metal- 
licity range will shed further light on the interpretation 
of these complex features. 

7.5. Fine Structure Emission Lines 

One of the puzzle s of the z ~ 3 LB G composite 
spectrum studied by IShaplev et al.1 (|2003| ) was the rel- 
atively strong emission from the excited fine structure 
lines Si II*, at 1265, 1309, and 1533 A. These lines are 
clearly visible in the spectrum of BX418 and in the z ~ 2 
composite; both spectra are shown in Figure[7l where the 
lines are marked. BX418 also shows significant emission 



from C II* A1335, which IShaplev et ail (|2003h did not 
detect, probably because of strong C II A1334 absorp- 
tion. Similarly, we see strong C II A1334 absorption in 
the z ~ 2 composite, but no C II* emission. Presumably 
the weakness of the low ionization interstellar absorp- 
tion lines in BX418 allows its detection here. These lines 
are not detect ed in the UV spectra of local starbursts, 
as discus sed bvlShaplev et al.| (|2003[ ) and shown more re- 
cently bv lSchwartz et al.ll2006l (see their Figure 13, which 
shows a clear comparison with the LBG composite spec- 
trum). Although the lines are weak, they are marginally 
stronger in BX418 than in the composite; this may be 
due to the stronger absorption in the composite, or it 
may suggest that the processes that give rise to these 
lines in high redshift but not local starbursts are even 
stronger in BX418. 
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iShaplev et al.l ((20031 attempted to reproduce these 
lines with Cloudy photoionization models but were un- 
able to do so, finding that all models which satisfacto- 
rily reproduced other observed lines underpredicted these 
lines by at least an order of magnitude. They concluded 
from this result that the lines are unlikely to originate in 
photoionized H II regions. We also attempt to use our 
Cloudy models to reproduce these lines, but similarly 
find that none of the models considered (whether they 
match the other properties of BX418 or not) come close 
to the observed line strengths. Even assuming that all 
of the emission from the Si II and C II multiplets arises 
at the excited fine structure wavelength, models which 
reproduce other observations underpredict the fine struc- 
ture emission by factors of 5-10 depending on the line. 
We conclude that if these lines do arise from the H II re- 
gions, either we are not modeling the regions correctly or 
some factor affecting these lines is not treated accurately 
by the models. 

The relative velocity of the lines may offer additional 
clues to their origin. In the z ~ 3 composite, these lines 
are redshifted by ~ 100 km s -1 , though they may be 
biased to the red by absorption on the blue side of the 
lines. The Si II* A1265, A1309 and C II* A1335 lines are 
also detected in the higher resolution spectrum of the 
lensed galaxy c B58, where they ar e again redshifted by 
~ 200 km s _1 (jPettini et al.l 12000ft . The lines are con- 
sistently redshifted in BX418 as well, by ~ 200 km s _1 , 
and in this case the weakness of the absorption makes 
significant biasing of the line centers less likely. This con- 
sistent redshift suggests an origin in the outflowing gas, 
perhaps through resonant scattering of photons originat- 
ing in the H II regions (though, as noted above, the pre- 
dicted emissio n from H II regions is not strong enough). 
Shapl ev" at al.l (|2003|) also considered the possibility of 
the lines arising in the outflow, but concluded that this 
explanation is not fully satisfactory because of the nar- 
rowness of the lines. The outflowing gas is apparently op- 
tically thin in these transitions, and therefore the lines 
should span the full range of velocities in the outflow, 
~ 1000 km s _1 or more. Instead they are quite nar- 
row. Similar arguments apply to BX418, where the Lya 
emission spans a range of over 2000 km s" 1 , but the fine 
structure lines are unresolved. However, it is also possi- 
ble that the line strength varies as a function of veloc- 
ity, and limited S /N prevents us from detecting emission 
except where it is strongest. The origin of these lines 
remains a puzzle. 

8. DISCUSSION 
8.1. Ionization Parameter 

A very high ionization parameter of log U ~ — 1 is sug- 
gestive of large numbers of massive, low metallicity stars 
forming in a small volume, in a configuration that must 
be considerably more extreme than anything observed 
locally. The source of the ionizing photons is unlikely to 
be an AGN, as indicated by the low C IV/C III] ratio and 
the nondetection of BX418 in deep 8 and 24 /xm images. 

The factors controlling the ionizati on parameter of a n 
H II region are discussed in detail bv lDopita et al.l (|2006f) . 
The ionization parameter increases significantly with de- 
creasing metallicity, because of the harder ionizing spec- 
tra and weaker winds of low metallicity stars. It also 



increases as the density of the H II region increases rela- 
tive to that of the surrounding medium, and has a weak 
inverse coupling with the pressure of the surrounding 
medium and a weak positive coupli ng with the mass of 
the star cluster. IDopita et al.l i|2()()(i: model the probabil- 
ity distributions of ionization parameters in local H II re- 
gions, finding a maximum value nearly an order of magni- 
tude l ower than the log U ~ —1 we infer for BX418. Sim- 
ilarly, iKewlev et al.l (|2001l) model the metallicities and 
ionization parameters of starbur st galaxies, and consider 
a max imum value of logt/ = — 2. Forstcr Schreiber et al.l 
(|2001[ ) have estimated the ionization parameter in the 
star- forming regions of the starburst galaxy M82, find- 
ing log[/ sa —2.3 on spatial scales of tens to hundreds of 
kpc. 

Although measurements of the emission lines needed 
to determine the ionization parameter are few at high 
redshift, the values inferred similarly do not approach 
that of BX418, although they are higher than the ioniza- 
tion parameters of l ocal galaxies of similar metallicities 
dPettini et all 120011 : iMaiolino et all 120081; iHainline et al.l 
I2009T ) . There is an exception, however, in the unusual 
emission line sp ectrum of the gravitationally lensed Lynx 
arc a t z = 3.36 (|Fosburv et al.H2003t IVillar-Martm et al.l 
l2004f h This object shows strong, narrow high ionization 
emission lines in the UV spectrum, including all of the 
lines observed in emission in BX418 and several others 
as well. Unlike BX418 it shows little to no continuum, 
so the interstellar absorption features we observe are not 
present. It is interpreted as a low mass H II galaxy un- 
dergoing a very young burst of low metallicity star for- 
mation, with an ionization parameter log U ~ — 1. Some 
of the line ratios are unlike those of BX418, however; in 
particular, C IV emission in the Lynx arc is 3.7 times 
that of H/3 (though this relies on a significant absorption 
correction), whereas we observe < 10% of H/3 at most. 
The spectrum of BX418 appears to be intermediate be- 
tween that of the Lynx arc and more typical high redshift 
galaxies. 

The difference in ionization parameter between BX418 
and that of more typical objects is more likely to be due 
to geometry than metallicity, since local galaxies with 
metallicities even lower than BX418 do not have ioniza- 
tion parameters nearly as high. We clearly have no data 
capable of constraining the geometry of star formation 
in BX418, so any conclusions we draw must be purely 
speculative. Given the scalings of the factors contribut- 
ing to the ionization parameter, however, the most likely 
scenario may be one in which the galaxy's massive stars 
are both numerous and densely packed enough that the 
galaxy is essentially a single large H II region with an 
extremely intense radiation field. The integrated spectra 
of more typical galaxies at z ~ 2 would then represent 
the superposition of individual H II regions, with proper- 
ties closer to those of the regions modeled locally. Clearly 
both additional data and additional modeling are needed 
to confirm or refute such a scenario. 

A potentially important implication of the extreme ion- 
ization parameter was raised in Section [721 The Cloudy 
photoionization models indicate that at very high ion- 
ization parameters and in the presence of dust grains, 
hydrogen-ionizing photons are selectively depleted by 
energy-specific grain opacities peaking near 912 A, re- 
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suiting in a decreasing flux in the hydrogen recombina- 
tion lines with both increasing ionization parameter and 
increasing metallicity, assuming that galax ies with higher 
meta llicities also have higher dust content (Bottorff et al. 
1998). Modeling suggests that the effect begins to be 
significant at log U ~ — 2, the minimum ionization pa- 
rameter required by our measurements of C III] and 
[O III]/[0 II] (these constraints do not depend on the 
decreasing flux of H/3, unlike that obtained from He II). 
The magnitude of the effect depends on both the amount 
and the size distribution of the dust grains, about which 
we know very little. We have not corrected the observed 
line ratios for extinction, given the low value of E(B — V); 
any correction applied would increase the fluxes of the 
UV lines relative to H/3, thereby requiring an ionization 
parameter even higher than that which we infer already 
We have, however, allowed grains to be included in the 
photoionization models (Cloudy includes Orion nebula 
grains as part of the H II region abundance set, and we 
have reduced them by varying factors to account for the 
low dust content of BX418). We have not attempted to 
use the models to constrain the dust properties of BX418, 
given the strong dependence on geometry, but we cannot 
reproduce the observed line ratios with entirely dust-free 
models. Such a situation is unlikely given the metallicity 
of ~ 1/6 Zq and the known correlation between metal- 
licity and dust content. 

For the grains included in the models, the hydrogen 
recombination line flux is reduced by a factor of ~ 2 
at log U ~ — 1. Such a reduction, while highly un- 
certain, would have significant implications. Most ob- 
viously, the true star formation rate would be a factor of 
~ 2 higher than the value we have inferred from Ha, and 
thus a factor of ~ 3 higher than the SFR obtained from 
the UV continuum. This is marginally inconsistent but 
not impossible; dust allowed by both the SED modeling 
and the Ha/H/3 ratio combined with a still rising UV 
continuum flux due to a young age could produce such 
observed SFRs, especially when combined with the sys- 
tematic and observational uncertainties. The generally 
good agreement of SFRs derived from Ha with those 
derived from all other av ailable diagnostics (jErb et al.l 
I2006ri iReddv et al.|[20Toh indicates that this effect can- 
not be widespread at z ~ 2, and indeed we do not expect 
it to be, given that typical z ~ 2 galaxies appear to have 
significantly lower ionization parameters. For BX418, 
more precise determinations of the star formation rate 
from a variety of indicators are needed to determine the 
magnitude of this effect. 

The reduction in Ha and H/3 flux would also affect any 
metallicity indicators employing those lines. In the case 
of the i?23 indicator used here (which takes ionization 
parameter into account) the Cloudy models suggest that 
the effect is mild, since the input metallicity of the best- 
fitting models is in reasonable agreement with the R23 
metallicity calculated from the predicted line ratios of 
those models. 

It is primarily the nebular He II emission which re- 
quires an ionization parameter as high as logf ~ — 1. 
This is uncertain; not only is the He II/H/3 ratio de- 
pendent on the diminution of the H/3 flux by grains, 
which are extremely difficult to constrain, the detec- 
tion of nebular emission relies on deconvolution with 
the broader stellar emission. As described in Section 



17.21 the two component fit we have adopted provides a 
better fit to the line at the 99% confidence level com- 
pared to a single, broad line, and it is obviously different 
in shape from purely stellar He II emis sion of similar 
strength (Dessaugcs-Zavadsky et al. 2010). Nevertheless 
we cannot rule out the possibility that the emission is 
purely stellar, with an unusual profile. However, even 
if the He II emission is stellar, there is significant ad- 
ditional evidence for an ionization parameter nearly as 
high: the lower limit on the [O III]/[0 II] ratio indicates 
> —2, as does the strength of the C III] emission, 
while the strong O III] emission indicates \ogU > —1.6, 
with only mild dependence on the diminution of H/3. We 
therefore consider the high ionization parameter to be 
a robust result, subject to the (not inconsiderable) un- 
certainties of the input ionizing spectrum used for the 
photoionization modeling. 

8.2. Comparisons with Other Galaxies at z ~ 2—3 

We have shown through comparison with the compos- 
ite spectrum of ~ 1000 z ~ 2 galaxies that BX418 is quite 
different from typical galaxies at this redshift. We have 
also constructed a composite spectrum of the youngest 
galaxies in the z ~ 2 sample for comparison; this spec- 
trum is generally quite similar to the composite of the 
full z ~ 2 sample, and the distinctive features of BX418 
are not appreciably stronger in the young subset of z ~ 2 
galaxies, suggesting that age alone is not enough to ac- 
count for BX418's uniqueness. A full study of the varia- 
tion of the spectra of z ~ 2 galaxies as a function of stel- 
lar population properties or particular spectral features 
is beyond the scope of this paper and will be presented 
elsewhere. However, it may be instructive to compare 
BX418 w ith the compos i te spe ctra of z ~ 3 Lyman break 
galaxies. Sha plev et al.l (|2003| ) divided the z ~ 3 sample 
into quartiles based on the strength of Lya emission, and 
found that, compared to the subset of galaxies with the 
strongest absorption, the strongest Lya emitters have 
bluer UV slopes, weaker low ionization absorption lines, 
slightly stronger He II emission, stronger emission from 
O III], and stronger emission in the C IV stellar wind 
line. These are all characteristics of BX418 (though the 
He II emission of BX418 is much stronger than that seen 
in any of the z ~ 3 composites), suggesting that BX418 
lies at the extreme end of a continuum of galaxy proper- 
ties that may be parameterized by the strength of Lya 
emission. 

Much has been written about the relationship of Lya 
emission to other galaxy properties and much remains to 
be understood, but perhaps the most reliable conclusion 
is that stronger L ya emission is associat e d with lower 
dust content (e.g. iVerhamme et al.l 120081 : iKornei et al.1 
120101 and references therei n), and perhaps mo re specif- 
ically with dust geometry. IKornei et al.l (2010) find that 
strong Lya emission is more common in older galaxies 
at z ~ 3, and propose that this is because supernova- 
driven outflows have cleared out channels through which 
the Lya photons can escape in the older galaxies. 

Other recent results are suggestive of the im portance 
of dust geometry as well. Red dv et al.l (|2010l ) have re- 
cently investigated the dust properties of galaxies at 
z ~ 2, using MIPS 24 /um imaging; they find that while 
galaxies with young best-fit ages have redder UV con- 
tinuum slopes, this does not necessarily mean that they 
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are dustier. Rather, they appear to have less dust than 
their UV slopes would imply, and are therefore better 
described by a steeper extinction law than more typical 
galaxies at z ~ 2 (BX418 has a bluer slope than most of 
these galaxies, and is inferred to have little dust content 
no matter which extinction law is applied). The rea- 
sons for this are not fully understood; younger galaxies 
may have a larger dust covering fraction, more consis- 
tent with a foreground screen, or the size distribution of 
dust grains may be different in younger objects. BX418 
is both young and blue, and we have inferred a low cov- 
ering fraction of neutral gas and dust; this may explain 
at least some of its differences from more typical young 
galaxies. While much remains to be clarified regarding 
the processes contributing to all of the features we have 
observed in BX418, it does seem clear that the covering 
fraction of gas and dust is a major factor in modulating 
the emergent UV spectrum of galaxies at high redshifts. 

8.3. Implications for Higher Redshifts 

We have suggested that because of its low mass, young 
age, low metallicity and low dust content, BX418 is likely 
to be similar to many of the galaxies now being discov- 
ered at 5 < z < 8. Detailed observations from future 
ground and space-based facilities are needed to confirm 
this claim, but if true the implications for high redshift 
galaxy formation and reionization may be significant. 
The Lya profile of BX418 shows substantial blueshifted 
emission, which is a signature of z ~ 3 galaxies with 
a significant Ly continuum escape fraction (Steidel et 
al., in preparation). If BX418 is indeed similar to typ- 
ical galaxies at higher redshifts, the implication is that 
the escape fraction may significantly increase with in- 
creasin g redshift, as alrea dy shown between z ~ 1 and 
z ~ 3 dSiana et al.l I2010T) and predicted at higher red- 
shifts (jWise fc Cenl l2009t iRazoumov fe Sommer-Larsenl 
[20TrT ). 

The detection of strong Lya and He II A 1640 emis- 
sion in a high redshift object is suggestive of the pres- 
ence of extremely low metallicity Population III stars, 
as discussed in de t ail by Uimenez fc Haimanl (|2006f) , and 
bv iFosburv et al. (2003) in the context of the Lynx arc. 
Given the many metal lines, our measured abundance of 
~ \ Z Q , and the fact that we are able to reproduce the 
observed line ratios with a normal stellar population, we 
do not believe this scenario applies to BX418. Indeed, 
models of Population III stars predict much higher He II 
equivalent widths of at least 5 < Wn Q n < 80 A, and Lya 
equiv alent widths of several hundred A as well ([Schaererl 
l200l . 

It is also worth noting that after Lya, the strongest 
emission line in the rest-frame UV spectrum of BX418 is 
not He II but C III] A1909. This line is observable from 
the ground to z ~ 10 and even beyond, and therefore 
future spectroscopic searches for galaxies at very high 
redshift may wish to target this line for unambiguous red- 
shift confirmation. Finally, study of BX418 may provide 
insight into the issue of correction of broadband magni- 
tudes to account for nebular line emission; we find that 
Ha emission contributes 0.4 ± 0.1 mag to the observed 
K-band magnitude, while combined [O III] AA4959,5007 
and H/3 emission account for 0.45 ± 0.4 mag in the H- 
band. Such corrections are important to consider when 



modeling the stellar populations of z ~ 7-8 galaxies, 
for which [O III] A A4959,5007 and K B fall in the Spitzer 
IRAC bands (e.g. lLabbe et al]|2010f ). Given the wider 
bandpasses of the IRAC bands, however, the correction 
for a similar galaxy at z ~ 7-8 would not be as large as 
those we find for the H and K bands. 

8.4. Remaining Questions and Further Work 

Perhaps the most important question we are unable 
to address with current observations is how many of the 
unusual spectral properties of BX418 are characteristic 
of galaxies with similar colors and metallicities, and how 
many are due to a very short-lived phase in its evolution. 
A larger sample of high quality observations of young, 
blue galaxies is needed to address this question, although 
this is a daunting prospect given their rarity and the 
substantial observing time required to obtain a spectrum 
of the S/N presented here. 

Local observations and improved modeling will also 
help to address this question. The expected strength 
of He II emission from W-R stars as a function of metal- 
licity and age is still not well known, and it is primarily 
the strong stellar He II emission that suggests that we 
may be observing BX418 at an unusual phase in its evo- 
lution. Similarly, the age and metallicity dependence of 
the P Cygni emission and absorption lines from massive 
stars is not well known; an improved understanding of 
these lines will be of great benefit to high redshift work, 
since they are among the stronger features in most spec- 
tra of distant galaxies. 

9. SUMMARY 

We have presented the analysis of the rest-frame UV 
and optical spectra of Q2343-BX418, a young (< 100 
Myr), low mass (M* ~ 10 9 M©), low metallicity (Z ~ 
1/6 Z ), and unreddened {E{B - V) ~ 0.02, UV contin- 
uum slope B — —2.1) galaxy at z = 2.3. In spite of its 
low mass, BX418 is an L* galaxy at z ~ 2, with a very 
low mass-to-light ratio (M/L)b = 0.03. Our primary 
conclusions are as follows: 

1. Using ratios of the strong optical emission lines, we 
measure an oxygen abundance of 12 + log(0/H) = 
7.9 ± 0.2, or Z ~ 1/6 Z©. We measure lower lim- 
its on the ratio of [O III]/[0 II], indicating a high 
ionization parameter, and of [N II] /Ha, indicating 
a low metallicity. We also determine the metal- 
licity via the direct, electron temperature method, 
using the ratio O III] AA1661, 1666/ [O III] A5007 
to determine the electron temperature and finding 
12 + log(0/H) = 7.8 ±0.1. 

2. The source of ionization in BX418 is unlikely to 
be an AGN. The C IV/C III] ratio is significantly 
lower than that seen in AGN, indicating a softer 
radiation field, and BX418 is undetected in deep 
imaging at 8 and 24 /^m, with a 24 /im upper limit 
well below the luminosities seen in most AGN of 
comparable redshifts and optical magnitudes. 

3. The rest-frame UV spectrum contains strong high 
ionization interstellar absorption lines from out- 
flowing gas, while the low ionization lines are ex- 
tremely weak. We compare the equivalent widths 
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and velocity profiles of the interstellar lines, and 
conclude that the covering fraction of high ioniza- 
tion gas is significantly higher than that of the low 
ionization gas, and thus that a larger than average 
fraction of the outflow is highly ionized. 

4. The Lya emission line (with rest-frame V^Lya = 54 
A) is extremely broad, with FWHM ~ 850 km s _1 , 
and shows significant blueshifted emission as well 
as the usual asymmetric redshifted profile. This is 
likely due to a low optical depth in the outflowing 
gas, which allows substantial escape of Lya pho- 
tons from the approaching gas on the near side of 
the galaxy. From the Lya/Ha ratio we find a Lya 
escape fraction f e = 0.4, assuming that Lya and 
Ha have the same spatial distribution. The true 
escape fraction may be higher than this, if Lya 
emission is more extended because of scattering. 

5. The UV spectrum contains strong He II A1640 
emission, which is well fit by the superposition of 
a broad stellar component and unresolved nebular 
emission. Photoionization modeling indicates that 
the nebular emission requires a very high ionization 
parameter \ogU ~ — 1. The strong stellar emission 
may require the presence of significant numbers of 
very young W-R stars. 

6. Strong nebular emission from C III] AA1907, 1909 
and O III] AA1661, 1666 is also present, and re- 
quires Z ~ 0.2 Z and \ogU > —1.6. We use the 
C III]/0 III] ratio to estimate the C/O abundance 
ratio, finding that it is consistent with the trend 
in C/O vs. O/H observed locally. Given BX418's 
young age and extreme star formation, this im- 
plies that the trend of increasing C/O with O/H 
is more likely to be due to metallicity-dependent 
stellar yields than to time delay effects, and that 
carbon production is not dominated by low mass 
stars. 

7. We observe strong P Cygni emission from the stel- 
lar wind lines C IV and Si IV, but are unable to 
reproduce both the emission and the absorption 
profiles with models. The strong emission requires 
a significant contribution from nebular emission, a 
stellar population age of only ~ 5 Myr, or a top- 
heavy initial mass function, while the absorption 
suggests a low metallicity of Z ~ 0.1 Z©. However, 
the similarity of the absorption profile to that of 
the higher metallicity composite spectrum of z ~ 2 
galaxies suggests that our understanding of this 
feature is incomplete. 

The properties of BX418 are broadly similar to those 
inferred for the young, low mass galaxies now being dis- 
covered at much higher redshifts, and therefore its de- 
tailed spectral properties may shed light on the likely 
physical conditions in these higher redshift objects, well 
before facilities are available to measure them directly. 
Important questions remain, however, that can be ad- 
dressed by observations both locally and at z ~ 2. A 
larger sample of high quality spectra of young, very blue 
galaxies at z ~ 2 is required in order to assess how 



many of BX418's unusual features are due to a short- 
lived evolutionary phase, and improved data and models 
of massive stars at low metallicity in the local universe 
are needed to aid in the interpretation of these spectra. 
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